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A  SPECTRAL  TRANSFER  MODEL  OF  OCEAN  WAVE  SPECTRA 


I.  INTRODUCTION 

In  this  repoit  we  develop  a  dynamical  description  of 
ocean  weve  spectra  that  can  be  applied  to  both  equilibrium  and 
developing  seas.  Our  principal  objective  is  to  develop  a 
model  that  can  be  applied  to  fairly  general  situations  so  that 
it  will  be  possible  to  correlate  data  taken  under  a  variety  of 
conditions . 

Although  the  most  casual  observation  of  ocean  waves 
shows  that  the  normal  ocean  is  cloeer  to  being  isotropic 
than  ont-dimensional,  we  will  develop  here  a  strictly  one¬ 
dimensional  model.  The  extension  to  two  dimensions  will  be 
treated  in  a  later  report. 

To  discuss  any  model  of  the  general  characteristics 
of  the  ocean  surface  as  we  do  in  the  following  sections ,  one 
must  have  a  clear  picture  in  mind  of  the  variety  of  wrvo 
structures  to  be  found  there.  For  the  sake  or  o3 arity  of  pre¬ 
sentation,  the  ocean  waves  are  separated  into  four  regions  of 
interest.  These  regions  are  determined  by  wavelength  and  have 
quite  different  characteristics.  The  longest  waves  are  the 
tides  and  tsunamis  ("tidal  waves")  which  have  a  wave  period 
of  from  10^  to  lO^  seconds  and  have  a  phase  velocity  on  the 
order  of  500  mph .  The  extremely  long  wavelengths  makes  these 


waves  uninteresting  in  the  present  di'  cussion  in  which  we 
are  concerned  with  rather  localized  wave  interactions .  The 
longest  waves  of  interest  to  us  are  called  swells  and  have 
periods  on  the  order  of  10  seconds  emd  wavelengths  between 
60  and  400  meters.  These  waves  are  sinusoidal  and  can  travel, 
without  loss  of  energy,  for  thousamds  of  miles  in  deep  water. 

It  is  through  these  long  wave  length  gravity  waves  that  the 
effect  of  storms  are  found  at  great  distances  from  storm 
centers . 

The  surface  gravity  waves  of  intermediate  size  ride  on 
these  much  longer  waves  and  can  interact  strongly  with  them. 

These  waves  are  seen  to  be  sharply  peaked  and  in  large  number . 
They  are  generated  by  the  wind  and  can  be  seen  to  grow,  in¬ 
crease  in  slope  emd  eventually  break  under  the  wind's  influence. 
The  generation  of  surface  waves  by  wind  is  an  incompletely 
understood  phenomenon,  but  the  mechanism  for  the  w’aves  breaking 
seems  to  be  clear.  A  wave  becomes  unstable  when  its  steepness, 
i.e.,  the  ratio  of  its  height  to  its  length,  exceeds  one  to  seven 
(Stokes  1847).  Therefore  a  wave  which  is  seven  meters  long 
will  break  when  its  height  exceeds  one  meter.  The  quantity 
H/X  ,  where  H  is  the  wave  height  and  X  the  wave  length, 
measures  how  close  a  given  wave  is  to  breaking.  A  quan¬ 
tity  we  will  use  frequently  is  the  wave  slope  (0)  defined 
as  0  *  ka,  where  k  is  the  wavenumber  and  a  the  wave 
amplitude. 
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The  smallest  waves  present  on  the  ocean  surface  are 
ripples  or  capillaries,  which  are  fractions  of  a  centimeter 
long.  These  waves  can  interact  strongly  with  the  surface  grav¬ 
ity  waves  in  the  neighboring  regime.  The  waves  in  this  region 
are  strongly  damped  by  molecular  viscosity  and  are  generated 
both  by  the  wind  and  by  the  gravity  waves. 

To  round  out  our  present  classification  of  the  ocean 
surface  effects,  we  must  introduce  the  notion  of  a  variable 
surface  current  which  can  arise  from  various  sources:  large 
scale  ocean  current  systems,  tides,  river  outlets,  wind-driven 
currents,  and  internal  wa  is.  An  internal  wave  is  unique  in 
that  it  gives  rise  to  a  travelling  current  pattern.  The  effect 
of  this  wave  at  the  ocean's  surface  is  to  increase  the  flow  of 
water  above  the  internal  wave  crest  and  decrease  the  flow  above 
the  internal  wave  trough,  thereby  giving  rise  to  a  position  de¬ 
pendent  translating  current  system  at  the  ocean's  surface.  Weak 
small  wavelength  surface  gravity  waves  can  interact  quite  strongly 
with  such  current  patterns  (see  Thomson  and  West,  1972) . 

To  capture  the  sequence  in  which  the  preceding  wave 
types  enter  a  region  of  the  ocean,  let  us  consider  a  glassy 
ocean  surface  in  which  there  are  no  currents  and  no  swells  from 
distant  regions.  If  a  uniform  wind  is  slowly  turned  on,  one 
observes  the  generation  of  ripples  on  the  water.  In  principle, 
it  is  the  shearing  stress  betv;een  the  wind  and  the  water  which 
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produces  this  initial  growth  of  the  wave's.  In  a  fairly  short 
time  interval,  these  ripples  grow  in  amplitude  and  become  in¬ 
termediate  size  gravity  waves.  This  time  is  dependent  on  the 
velocity  of  the  wind.  Waves  of  many  sizes  are  generated;  those 
with  sma’l  wavelength  are,  however,  quickly  damped  by  viscosity. 
Those  with  longer  wavelength,  as  they  grow,  become  more  effi¬ 
cient  in  extracting  energy  from  the  wind  since  their  amplitude 
acts  as  a  sail  on  which  the  wind  can  push.  The  waves  which 
are  the  most  efficient  in  extracting  energy  from  the  wind  are 
those  which  have  a  phase  comparable  to  the  wind  speed. 

The  continued  growth  of  these  primary  gravity  waves 
is  of  course  limited  by  the  stedaility  criteria  mention  earlier, 
i.e.,  H/X  »  1/7.  The  energy  being  supplied  by  the  wind  must, 
therefore,  be  redistributed  to  other  wavenumbers.  The  inter¬ 
action  between  gravity  waves  is  weak,  however,  so  energy  can 
be  fed  in  by  the  wind  more  rapidly  than  it  can  be  redistributed 
by  the  interaction  between  waves.  This  excess  energy  is 
dissipated  in  two  ways;  (i)  the  formation  of  "white  caps"  in 
which  the  gravity  wave  becomes  unstable  and  breaks  and  (ii)  the 
generation  of  "parasitic  capillaries",  i.e.,  short  waves  which 
are  radiated  from  the  crest  of  a  sharply-peaked  gravity  wave 
emd  are  subsequently  viscously  damped. 
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'  If  the  wind  speed  is  insufficient  to  cause  the  sharp 
cresting  of  the  gravity  waves  which  lead  to  (i)  and  (ii)  above, 
these  effects  may  still  be  induced  by  means  of  a  number  of 
diffe::ent  mechanisms.  The  "capping"  or  breaking  of  waves 
can  be  caused  by  two  fully-developed  waves  running  together 
so  that  their  superposition  forms  a  sharply-crested  peak  which 
is  unstable.  Breaking  can  also  occur  by  relatively  short  gravity 
waves  passing  over  the  crest  of  a  swell  so  that,  again,  the 
crest  of  the  combined  waves  is  unstable.  It  should  be  pointed 
out  that  breadcing  and  the  formation  of  "parasitic  capillaries 
are  competitive  mechanisms  and  where  one  is  found  so  generally 
is. the  other.  It  is  probably  the  case,  however,  that  the  latter 
mechanism  is  operative  when  one  exceeds  the  stability  limit  by 
only  a  small  amount. 

As  one  increases  the  wind  speed,  longer  wavelength  waves 
are  generated.  Those  waves  already  present  become  saturated 
and  begin  breaking.  The  energy  pumped  into  the  waves  is  rapidly 
dissipated  by  "capping"  and  the  fraction  of  those  waves  which 

are  breaking  increases  rapidly. 

There  have  been  a  number  of  attempts  to  develop  a  sta¬ 
tistical  description  of  the  air-sea  interface.  The  most  ex¬ 
tensive  modeling  of  this  type  has  been  carried  out  in  a  series 
of  papers  by  Hasselman  (1962,  1963). 
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The  basic  concept  behind  Hasselman's  approach  stems 
from  the  fact  that,  although  the  equations  describing  the 
wave  motion  at  the  sea  surface  are  non-linear,  these  inter¬ 
actions  are  remarkably  weak  except  under  extreme  circumstances . 
Perturbation  theory  is  then  appliceible.  Under  these  conditions 
transfer  of  energy  occurs  most  rapidly  between  various  modes 
that  satisfy  a  resonance  condition,  fiasselman's  model  has  been 
applied  successfully  to  the  description  of  wave  spectra  at 
relatively  long  wavelengths  where  the^  very  strong  non-linearities 
that  result  in  breaking  waves  zuid  white  caps  have  only  an 
indirect  effect. 

However,  in  the  normal  sea,  a  significant  fraction  of  the 
wave  spectrum  is  usually  assumed  to  be  "saturated".  In  this 
region  the  continued  energy  input  from  the  wind  is  balanced  by 
highly  non-linear  effects  eventually  resulting  in  dissipation 
of  the  energy.  These  brezdcing  waves  or  white  caps  are  not 
described  by  the  Hasselman  model  (or  any  model  based  on  per¬ 
turbation  theory . ) 

The  approach  we  use  here  is  wO  construct  a  "trrmsfer" 
equation  in  which  different  interaction  effects  are  identified 
emd  models  constructed  for  each  of  these  effects.  This  transfer 
equation  and  tJie  derivation  of  the  individual  terms  are  des¬ 
cribed  in  detail  in  the  next  section. 
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2.  THE  TRANSFER  EQUATION 

We  wish  to  develop  a  phenomenological  equation  which 
describes  how  energy  is  fed  into  ocean  waves,  how  it  gets  re¬ 
distributed  among  these  waves  and  how  it  is  dissipated.  To 
realize  this  goal,  it  is  necessary  to  summarize  in  a  single 
equation  all  the  effects  discussed  qualitatively  in  the 
preceding  section.  In  general,  we  need  to  construct  an  equation 
for  the  rate  of  change  of  wave  energy  at  a  given  poiot  in 
space  (x) ,  time  (t) ,  and  wavenumber  (k) .  Since  we  are  particu¬ 
larly  interested  in  constructing  a  model  that  allows  for  spatial 
gradients  resulting  from  a  travelling  current  field,  we  think 
of  a  transfer  equation  that  will  allow  waves  to  be  born  at 
some  point,  to  convect  to  a  different  point  in  space,  and  to 
decay  or  dissipate.  To  this  end  we  have  divided  the  present 
section  into  a  number  of  subsections,  each  one  dealing  with  a 
single  term  in  our  equation.  Figure  1  gives  a  pictoral  summary 
of  the  effects  we  include. 

In  Figure  1  we  divide  the  waves  on  the  ocean  surface 
into  three  regimes  as  mentioned  previously.  The  arrows  in  the 
figure  indicate  the  direction  of  energy  flow  between  different 
wavenumbers;  the  ledaels  on  the  arrows  indicate  the  mechanism  by 
which  the  energy  is  tremsferred.  The  solid  arrows  indicate 
that  the  mechanism  is  important  in  the  wavenumber  region  con¬ 
sidered;  a  dashed  arrow  indicates  a  relatively  unimportant 
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mechanism.  For  example,  viscosity  is  an  important  dissipative 
mechanism  for  high  wavenumber  (k  }  waves  but  not  for  low  wave- 

n 

number  waves  (k  )  . 

Lj 

We  define  the  quantity  F(k,x,t),  [sometimes  written  simp 
as  F()c)  for  brevity]  ,  as  the  spectral  density  of  the  wave  ampli 
tude  at  the  position  x  at  time  t  for  waves  with  wavenumber  k. 
This  is  the  quantity  that  describes  the  processes  represented  i 
Figure  1.  The  total  rate  of  change  in  F(k,x,t)  is  dependent  on 
the  five  processes  labeled  in  Figure  1:  (1)  wind,  (2)  resonant 

interactions  between  gravity  waves,  (3)  the  radiation  stress  of 
a  wave  from  one  k-regime  passing  over  a  wave  of  another  regime, 
(4)  viscous  damping,  (5)  formation  of  parasitic  capillaries 
and  wave  breaking.  The  overall  form  of  the  one-dimensional 
equation  which  we  will  arrive  at  is  shown  below: 


3F 


=  -fc  +tn  11.  -  ^  H 

3t  g  ^  3x  3U  dx  3k 


r  ^  F  +  1-  /n 
dx  3k  \  k  Tic  j 


(»).„•«) 


•  (s) 

breaking  '  'dissipation 


The  first  four  terms  on  the  right  hand  side  comprise,  respec¬ 
tively,  the  effects  of  convection,  refraction,  radiation  stress 
and  wave-wave  interaction.  In  the  remainder  of  this  section 
we  discuss  each  of  the  terms  of  this  equation  in  detail. 


2 . 1  Model  of  the  Total  Time  Derivative 

We  will  restrict  our  attention  in  the  following  discus¬ 
sions  to  tha  one-dimensional  case.  We  write  the  energy  den¬ 
sity  at  a  point  on  our  one-dimensional  ocean  surface,  where 
there  exists  a  current  [U(x)],  as  a  function  of  the  wavenumber 
(k)  ,  position  (x)  ,  and  time  (t)  .  To  c’etermine  the  total  rate 
of  change  in  the  energy  density  at  a  point  in  space  one  must 
consider  how  the  energy  density  changes  v;ith  respect  to  all  its 
arguments  in  an  increment  of  time  At. 

The  quantity  F(k,x,t)  is  the  probable  value  of  the 
power  spectrum  of  the  wave  amplitude  at  the  phase  point  k,x,t 
and  is  proportional  to  the  energy  density.  Its  rate  of 
change  with  time  dF/dt  combines  several  effects.  In  the 
following,  we  will  derive  models  for  the  rate  at  which  this 
energy  density  will  change  as  viewed  by  an  observer  who  is 
being  convected  along  with  the  individual  wave  packets.  The 
trajectory  of  a  packet  in  our  two-dimensional  phase  space  is 
specified  in  terms  of  the  local  dispersion  relation  u)(k,x) 

(we  only  consider  here  the  special  case  of  a  time  independent 
dispersion  relation) .  The  generalization  to  the  time  depend¬ 
ent  case  [w»a)(k,x,t)]  will  be  made  subsequently.  The  trajec¬ 
tory  along  which  energy  naturally  prcpagates  is  then  given  by 


and 


dt 


dk 

It 


V..  0)  or  in  one  dimension  srr 
k 


dO) 

9k 


-  oj  or  in  one  dimension 

X 


dk 

dt 


9a) 

9x 


(2.1.1) 


When  0)  is  explicitly  a  function  of  a  local  surface  current 
[U(x)]  this  latter  equation  becomes 


dk 

dt 


9U  dx 


in  one  dimension. 


(2.1.2) 


For  the  sake  of  maintaining  a  consistent  convention  with  other 
writers,  we  express  dx/dt  in  the  form 


dt 


C  +  U 

g 


(2.1.3) 


where  is  the  group  velocity  relative  to  the  local  water  sur¬ 
face  (Cg  =  9a)^/9k  where,  for  gravity  waves,  0)^  =  /^).  Also 
for  a  steady  current 


a!  =  Uk  +  u^(k)  =  constant 

and  using  Equation  (2.1.2)  we  have 

dk  _  .  dU 
dt  dx 

and  from  Equation  (2.1.4), 

k(C+U)  =  0)  =  constant 


(2.1.4) 


(2.1.5) 


where  C  is  the  phase  velocity  in  the  ocean  coordinate  system 
(C  =a)^/k)  . 


11 


The  total  rate  of  change  of  F(k,x,t)  following  individual 
wave  packets  is  related  to  the  local  space- time  gradients 
according  to 

dt  3t  dt  3x  dt  ck  gain  ^loss 


or,  using  Equations  (2.1.3)  and  (2.1.5), 


(2.1.6) 


^  =  31+  (c  +u)^-k^^ 

dt  3t  '  g  3x  dx  3k 


F  .  -  F, 

gain  loss 


The  first  two  terms  on  the  left-hand  side  of  Equation  (2.1.6) 
are  familiar,  being  the  variation  in  F  with  time  and  that  due 
to  convection  at  a  point  in  physical  space.  The  third  term  is 
a  little  more  interesting,  since  it  results  from  the  variation 
in  wave  number  due  to  the  gradient  of  the  ocean  current  (refrac¬ 
tion)  .  This  term  is  analogous  to  the  action  of  an  external 
force  on  a  group  of  particles  undergoing  collisions  (F  would  be 
the  number  of  particles).  Equation  (2.1.6)  is  the  wave  analog 
of  the  Boltzmann  equation. 

In  problems  connected  to  the  surface  currents  produced 
by  propagating  internal  waves,  the  coordinate  system  in  which 
the  wave  pattern  is  steady  is  that  moving  with  the  internal  wave 
at  its  phase  velocity  (Cj)  .  It  will  be  convenient  in  these 
cases  to  write  the  total  current  seen  in  this  coordinate  system 
in  the  form 

U  =  AU  -  Cj 

where  AU  is  the  flow  perturbation  produced  by  the  internal  wave. 


12 


2.2  Model  for  VJave  Breaking 

As  discussed  previously,  there  are  three  types  of 
non-resonant  interactions  that  can  give  rise  to  wave  breaking: 
(i)  the  superposition  of  two  fully-developed  waves,  (i:.)  the 
passage  of  a  packet  of  short  waves  over  a  fully-developed  wave 
of  low  wavenumber  and  (iii) ,  the  continued  action  of  the  wind 
on  a  wave.  The  first  two  of  these  interactions  occurs  rapidly, 
typically  on  a  time  scale  of  the  order  of  the  wave  period  or 
less.  The  action  of  the  wind  occurs  over  a  longer  time  period. 
In  this  section  we  develop  a  model  for  the  breaking  of  fully- 
developed  waves. 

The  basic  model  we  will  use  is  dependent  on  there  being 
some  critical  value  for  the  ratio  of  wave  height  to  wave 
length.  The  analysis  of  Stokes  established  that  for  gravity 
waves  this  value  was  approximately  1/7  after  which  gravita¬ 
tional  instability  would  set  in.  Referring  to  this  ratio  as 
the  wave  slope,  we  assume  that  if  the  wave  slope  exceeds 
some  critical  value  (©q) r  strong  interactions  occur  which 
radically  alter  the  energy  distribution  in  wavenumber  space. 
These  mechanisms  could  be  wave  breaking  or  the  generation  of 
parasitic  capillaries  at  the  wave  crest.  Exactly  how  these 
effects  occur  is  irrelevant  here;  we  simply  assume  that  when¬ 
ever  the  slope  (0)  attempts  to  exceed  0^,  some  strongly 
dissipative  process  occurs  to  keep  0  at  0^.  The  problem  which 
remains  is  to  estimate  the  proper  statistics  of  the  wave  slope. 


We  assume  that  the  wave  field  of  the  ocean  is  a  random 
superposition  of  modes;  further,  that  thei^e  waves  have  random 
phases  and  are  distributed  over  a  broad  range  of  wavenumbers. 
If  we  write  the  slope  of  the  surface  0  as  the  sum  of  the 
contributions  from  each  mode,  then 


0 


(2.2.1) 


By  assumption  each  of  the  0j^  is  a  random  variable,  so  that  if 
no  non-linear  process  existed  in  the  interaction  between  the 
surface  waves,  the  law  of  large  numbers  would  result  in  a 
Gaussian  distribution  for  the  surface  slope.  In  that  case  it 
would  be  possible  to  write  for  the  probability  that  the  sur¬ 
face  has  a  slope  greater  than  the  critical  value  (0^)  as 


P(|0|  >  0^)  =  de  exp  [-  0V^  }  . 


(2.2.2) 


7T0^  '<> 


In  the  asymptotic  region  j^0^  >> 
approximate  expression 


^0^  >>  ^  j ,  we  may  write  the 


(2.2.3) 


to  replace  Eg.  (2.2.2).  Eq.  (2.2.3)  is  the  fraction  of  the 

ocean  surface  in  the  wave  field  whose  slope  exceeds  the  critical 

value  0  and  would  form  white  caps  according  to  this  model, 
o 


y* 
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The  coherence  'cime  or  duration  of  the  breaking  pro¬ 
cess  is  of  the  order  of  the  recifiocal  of  the  wave  period, 
i.e.,  the  wave  formed  by  the  superposition  of  the  surface 
modes.  If  a  fraction  e  of  the  wave  energy  is  lost  in  one 
such  breaking,  then  the  mean  rate  of  loss  of  the  energy 
density  is  of  the  order 


(  =  -(*)£  P(  |0l  >  0^)  F  .  (2.2. 

breaking  ° 

Equation  (2.2.4)  is  the  probable  fraction  of  ej-iergy  lost  per 
unit  time  interval  due  to  wave  breaking.  Using  (2.2.3),  we 
may  rewrite  (2.2.4)  as 


(H) 


=  -u)e  F 


breaking 


0. 


{</^^] 


(2.2.5) 


in  terms  of  the  mean  square  slope  of  the  surface  0  . 

The  mean  square  slope  0  may  be  expressed  in  terms  of 
the  energy  density  by  recalling  that 


/: 


F(k)dk  =  a^ 


(2.2.6) 


where  is  the  mean  square  elevation  of  the  ocean  surface.  We 
may  use  this  to  write  the  mean  square  slope  as. 


- T  f° 

=  (VC) '  =y 


k'^  F(k)dk  . 


(2.2.7) 


1- 
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In  order  to  develop  the  concept  of  a  spectral  model  for  non- 
linearly  interacting  waves,  it  is  necessary  to  invoke  the  con¬ 
cept  of  a  correlation  interval  in  wavenumber  space.  We  assume 
that  waves  within  the  interval  k  ±  Ak/2  arc  completely  (phase) 
correlated  and  act  as  a  single  non-linear  wave.  Waves  separ¬ 
ated  by  a  wavenumber  interval  greater  than  Ak  are  assumed  to 
be  uncorrelated  with  no  consistent  phase  relationships.  Since 
for  a  well-developed  sea  we  expect  Ak  to  be  proportional  to 
k,  we  take  the  mean  square  slope  of  waves  havxng  wavenumber 
k  to  be  _ 

0^(k)  =  constant  •  k^F(k)  .  (2.2.8) 

2 

The  constant  can  be  absorbed  into  the  factor  0  and  we  can 

o 

write  the  order  of  magnitude  expression. 


We  may  extend  this  result  to  the  case  where  we  treat 
the  right  and  left  traveling  waves  separately.  This  separation 
occurs  when  we  view  the  oceam  in  a  coordinate  system  in  which 
the  local  current  is  steady.  In  this  coordinate  system,  the 
probability  of  the  net  slope  exceeding  0^  in  the  wave  should 
depend  on  the  total  energy  density  F  =  +  F^  ^ ,  Thus, 

we  have  the  set  of  coupled  equations. 
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and 


F 

3t 

breaking 

J 

breaking 


(2.2.10) 

(2.2.11) 


2 . 3  Non-Linear  Wave-Wave  Interaction 

In  one  sense  the  non-linear  interaction  between  sur¬ 
face  gravity  waves  is  the  most  clearly  defined  probleni 
discussed  in  this  report  in  that  there  exists  a  well-defined 
system  of  equations  which  describes  the  interaction  process. 
These  equations  are  the  kinematic  boundary  condition  at  the 
ocean  surface  and  Bernoullis'  equation.  Solving  these 
equations  for  a  surface  profile,  for  any  but  the  simplest 
situation, is  a  notoriously  difficult  task.  One  must,  there¬ 
fore  turn  to  suitable  approximation  techniques  to  solve 
these  equations.  The  approach  which  is  the  most  compatible 
with  the  spectral  modeling  developed  in  this  report  is  a 
modal  description  of  the  interaction  process. 

The  mode  oriented  models  treat  the  ocean  surface  as 
a  superposition  of  waves  (in  the  unperturbed  ocean  this 
would  be  a  superpositon  of  sine  waves).  The  non-linear  terms 
in  the  interaction  equations  for  the  surface  waves  will 
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allow  coupling  of  the  different  modes.  The  first  systematic 
development  of  such  a  theory  is  that  due  to  Hasselmann  (1962) 
and  (1966) .  A  similar  approach  has  been  used  by  West,  Watson 
and  Thomson  (1973)  to  construct  an  interaction  equation 
in  terms  of  the  slope  of  a  wave  qj^  , 


dq 


(s) 


dt 


E 


ki+k2-=k 


(Sj^)  (S2)  i  (Sj^a)^+S2W2-SuJj^)  t 

^1,  e 


k,k2  ^kj^  ^k2 


kk3  (Sj^)  (S2)  (S3)  ^“k^k^kk  ^ 

^k  k  ^k  ^k  ^k  ® 
k,+k.=k.+k  ''l''2  ‘'l  ^^2  ^3 

^  ^  (2.3.1) 


where 


^k^k2kk3  ~  ^l“kj^  ^2“k2  ^3‘^k3  "  ®“k 


(2.3.2) 


and  s=±  corresponds  to  a  right  (+)  or  left  (-)  travelling 
wave  and 


=  +1  if  s^  *  + 


in  the  exponentials.  In  the  language  of  Hasselmann  (1966)  , 
Eq.  (2.3.1)  includes  three  and  four  wave  scatterings  in 
a  perturbation  expansion  of  the  wave  amplitudes. 

The  structure  of  the  interaction  terms  in  Eq.  (2.3.1) 
indicate  that  strong  interactions  and  energy  transfers  occur 
between  modes  which  together  form  a  resonance,  i.e.. 


0) 


~ 


ki  +  k2  =  k 


2nd  order 


or 


(ill  +  (1)2  =  ^2  ■*■ 


ki  +  k2  =  ks  +  k 


1 

/ 


3rd  order 


For  gravity  waves  (oj^  =  /gk^)  ,  it  is  well  known  that  the  only 
resonant  solutions  at  2nd  order  are  the  uninteresting  ones 
k^  or  k2  =  0.  At  third  order,  however,  non-zero  solutions 
are  possible.  In  order  to  have  a  pertect  frequency  match, 
the  only  solutions  permissible  are 
ki  =  k^  and  k2  =  k 


or 


k^  =  k  and  k2  =  k^ 


where  k  (=  k^)  is  the  fourth  member  of  this  interacting  set  of 
waves.  The  near  resonance  condition  is  depicted  in  Figure  (2) 
where  a  four  wave  interaction  is  represented  schematically  in 
(2a)  and  the  location  of  the  waves  on  the  linear  dispersion 
curve  is  indicated  in  (2b)  . 

We  see  in  Figure  (2b)  that  a  slight  detuning  has  been  allowed 
to  occur.  Fortunately,  for  finite  energy  transfer  times  it  is 
not  necessary  to  demand  a  perfect  frequency  match,  and  modes 
which  couple  together  to  give  a  frequency  mismatch  601  <  1/T 
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(where  T  is  the  effective  lifetime  for  the  energy  transfer) 
will  form  effectively  a  resonant  group.  The  magnitude  of  the 
mismatch  frequency  may  be  expressed  in  terms  of  the  devia¬ 


tion  of  the  wavenumbers  and  <2  from  our  reference  wavenum¬ 
ber  k : 

=  k(l  +  6^) 


k2  =  k(l  +  62) 


(2.3.3) 


The  wavenumber  k^  is  a  function  of  kj^  and  k2*. 

k^  =  kj^  +  k2  -  k  .  (2.3.4) 

Then,  to  second  order  in  6^^  or  have  for  the  mismatch 

6u)  =  +  U)2~W2"'W 


The  region  of  wavenumber  space  that  contributes  to  the  resonance 
(600  <  1/T)  is  sketched  in  Figure  3  for  wT  ~  250  (this  value  of  uT 
is  of  the  order  expected  in  a  near  saturated  sea) .  Thus,  for 
the  small  values  of  Sm  that  contribute  to  the  integral,  the 
strongly-interacting  modes  _ie  close  together.  This  fact  can 
be  used  to  derive  a  relatively  simple  description  of  the  wave- 
wave  interaction.  To  show  this  we  start  with  the  non-linear 
mode  interaction  formalism  of  Hasselman  (1962) .  He  uses  per¬ 
turbation  theory  to  describe  the  resonant  energy  transf'^r  be- 
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FIGURE  3. 


The  Region  of  Wavenumber  Space  that 
Contributes  to  the  Resonance  (ow  <  1/T) 
is  Sketched  for  wT  ~  250 


tween  different  inodes  of  the  ocean  spectrum.  By  assuming  the 
statistical  distribution  of  wave  amplitudes  on  the  ocean  sur¬ 
face  tc  be  Gaussian,  he  was  able  to  derive  an  expression  for 
the  rate  of  increase  of  energy  of  a  particular  mode  (k^)  in 
terms  of  the  product  of  the  energy  density  in  three  other 
modes 


^  =  A(“4"i 


F2F3  +  -  a)2F 


lF,F4-.jF2F3F^)d\d2k. 


where 


(2.3.6) 


2  2 

9irg^  D  w. 

A  =  — - — - ^ - —  6  (u),  +  W-  -  oj,  -  u) . ) 

A  1  \2  12  3  4 

4p  (0)^0)20)300^) 


and  we  have  required  that 


In  Eg.  (2.3.6)  F^  represents  the  power  spectrum  of  the 
wave  height  at  the  wavenumber  [F^  =  F(ic^)].  According  to 
our  earlier  discussion,  the  major  contribution  to  the  integral 
will  come  from  wavenumber  and  k2  that  are  close  to  k^ 

(see  Figure  3).  If  this  interval  is  small  enough,  we 
are  able  to  expand  F(k^)  in  a  Taylor  series  about  k^  =  k^  and 
express  the  integral  in  terms  of  the  derivatives  of  F(ic)  at 
k  =  .  To  do  this  explicitly,  we  would  need  to  proceed  beyond 

the  limits  of  perturbation  theory  to  derive  specific  procedures 
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for  allowing  for  the  frequency  mismatch.  However,  the  form 
of  the  result  can  be  obtained  without  requiring  knowledge 
of  such  methods. 


In  general,  we  will  have 


at 


G 


3F  a^F 

W  '  ••• 


(2.3.7) 


We  include  only  terms  up  to  the  second  derivative  and  note  that 
this  interaction  can  only  rearrange  the  energy  distribution, 
leaving  the  total  energy  constant .  Thus 


(2.3.8) 


Since  this  relation  must  hold  for  arbitrary  distributions, 
3F 

F(k) ,  ^  must  be  a  derivative  of  the  form 


aF 


a 


|^g(k,F, 


(2.3.9) 


Using  the  general  structure  of  Eq.  (2.3.6)  we  have  that  the 
right-hand  side  of  Eq.  (2.3.9)  must  be  cubic  in  F,  so  that 
Eq.  (2.3.9)  must  be  of  the  general  form 


(constant)  j  . 


(2.3.10) 


Dimensional  analysis  can  be  used  to  further  limit  the  parameters 
for  gravity  waves.  Since  the  acceleration  of  gravity  is  the  only 
available  parameter  and  since  F  has  dimensions  of  (length)^,  we 
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can  only  have 


ar 

at 


Cj^Ti  /g 


ak 


,  8.5-s 
k 


ak 


(Fk®) 


(2.3.11) 


where  is  a  dimensionless  number. 

Finally  we  argue  that,  since  there  is  no  characteristic 
scale  to  the  problem,  an  equilibrium  condition  must  correspond 
to  one  in  which  the  wave  slopes  are  independent  of  scale, 
i .e . ,  that 


(Fk^) 


equilibrium 


constant 


Thus,  we  choose  s  =  3  and  write 


(2.3.12) 


wave -wave 


a 


(k^F)^ 


(2.3.13) 


The  value  of  can  only  be  determined  from  a  more  de¬ 
tailed  theory  or  from  experiment,  such  as  to  that  carried  out 
by  Benjamin  and  Feir  (1967).  Here  a  finite  amplitude  travelinc 
wave  was  set  up  in  a  wave  tank  and  its  decay  or  break  up  observed 
as  a  function  of  distance  from  the  wave  maker .  In  the  transfer 
model,  the  wave  height  should  obey  in  this  case 

f  Vi" 

*  The  Ben jconin-Feir  study  consisted  of  the  observation  of  the 
breakup  of  a  single  coherent  non-linear  wave.  The  appro¬ 
priate  experiment  to  test  a  diffusion  model  should  be  some¬ 
what  different  and  would  observe  the  decay  of  a  coherent 
wave  in  the  presence  of  an  already  fully  developed  sea. 


It  is  convenient  for  calculational  purposes  to  define  a 

3 

(dimensionless)  quantity  y  »  Fk 


(2.3.15) 


This  equation  can  be  interpreted  as  a  diffusion  equation  (in 
wavenumber  space)  having  a  non-linear  diffusion  coefficient 
defined  by 


2 

K  Y 


(2.3.16) 


An  effective  measure  of  the  time  (T)  for  energy  to  spread 
out  of  a  wavenumber  interval  Ak  (i.e.,  a  lifetime)  is 


(2.3.17) 


We  argue  that  the  permissible  mismatch  of  frequency  6a)  must 
be  of  order  T  ^ .  Thus  we  may  write 


6a) 

0) 


(k3-k^) 


a)T 


4  k' 


‘^k 
■"  ■ 

a)Ak^ 


(2.3.18) 


Since  both  k^^-k^  and  k^-k^  are  each  of  order  Ak^, 
we  find  that  the  width  of  the  wavenumber  regime  into  which 
the  mode  decays  is  of  order 


(2.3.19) 


I 

) 

i 

t 

it 

tl 
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The  lifetime  is  giver*  by 


wT  =  - = — 

^■^T\  y 


(2.3.20) 


For  a  saturated  sea,  the  value  of  y(= 
experiment : 


is  known  from 


Y  .  a  0.00585  ±  10% 

S  d  u 


(2.3.21) 


Thus  tb^  lifetime  (N,  measured  in  wave  periods)  of  a  mono¬ 
chromatic  wave  in  a  near  saturated  ocean  is  expected  to  be 
given  by 


wT 


(.00585) 


30.5 


(2.3.22) 


Thus,  experiments  of  this  type  could  be  used  to  determine  the 
coefficient  Cj^,  as  well  as  simply  to  test  the  diffusion  .lodel. 


2. 4  Linear  Interaction  Theory 

2.4a  Model  of  the  Surface  Wave-Current  Interaction 

In  this  section  we  will  be  concerned  with  constructing 
a  model  to  describe  the  effect  of  long  waves  on  short  waves; 
i.e.,  the  interaction  between  widely  separated  regions  of 
wavenumber  space.  These  interactions  may  be  separated  into 
a  number  of  distinct  physical  situations,  such  as  swells 
passing  through  a  region  of  developing  ocean,  surface  currents 
produced  by  internal  waves  interacting  with  surface  waves,  etc. 

To  establish  a  model  of  these  phenomena,  we  will  review  some 
aspects  of  the  linear  theory  for  these  processes.  We  consider 
the  short  waves  to  be  infinitesimal  surface  waves  in  the 
development  of  the  linear  theory  and  the  long  waves  to  be 
describable  in  terms  of  a  surface  current.  The  analysis  centers 
on  how  energy  is  transferred  between  the  two  wave  length  regimes. 
The  model  is  an  outgrowth  of  a  series  of  studies  made  by  Longuet- 
Higgins  and  Stewart  (1961)  and  (1964)  and  more  recently  by 
Phillips  (1966,  pages  50  and  136). 

The  model  ocean  used  in  this  section  is  both  homogeneous 
and  irrotational.  This  allows  us  to  use  a  potential  description 
of  the  velocity  field,  i.e. 

V  X  u  *  0  which  implies  u  =  V*  .  (2.4.1) 

The  velocity  potential  also  satisfies  Laplace's  equation 

V^4>  =  0  (2.4.2) 
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since  the  fluid  is  assumed  to  be  incompressible  (V*u  =  0). 

Using  the  integral  of  the  momentum  equation,  we  may  write 
Bernoulli's  equation  at  the  ocean  surface  as 

1  +  gz  =  0  at  z  =  C  (2.4.3) 

with  the  condition  that  hydrostatic  equilibrium  must  prevail 
at  infinity,  and  z  =  C  is  the  free  surface  of  the  ocean.  The 
rate  of  increase  in  the  wave  height  following  a  fluid  element 
is  the  vertical  component  of  the  fluid  velocity: 

^  .  1^.  at  J  -  C 

dt  3z 

where 

^  =  I—  +  V<I)*V  (2.4.5) 

3t  at 

is  the  Eulerian  derivative. 

We  intend  to  deal  with  the  linearized  forms  of  Equations 
(2.4.3)  and  (2.4.4)  and  so  we  separate  the  velocity  potential 
and  surface  displacement  into  a  long  and  short  wave  component; 

+  (2.4.6) 

and 

5  =  X  +  h  . 

The  terms  t  and  x  satisfy  the  free  surface  equations  independently 
of  the  surface  wave  height  (h)  and  velocity  potential  (((>), 


that  is 


^t  +  j  VV  •  VV  +  gX  -  0  at  z  -  X  (2.4.7a) 

and 

+  VV-VX  =11  at  z  =  X  .  (2.4.7b) 

These  equations  refer  to  the  unperturbed  ocean  surface;  that  is, 
unperturbed  by  the  short  surface  waves  and  may  be  identified 
with  swells,  surface  currents  produced  by  internal  waves  or  the 
effect  of  other  generating  mechanisms  of  interest. 

Equations  (2.4.7ci)  and  (2.4.7b)  are  valid  at  the  unperturbed 
surface  z  ■  x»  but  the  more  general  expression  given  by  Equations 
(2,4.3)  and  (2.4.4)  are  valid  at  the  free  surface  z  ■  X  +  h.  To 
determine  the  interaction  between  the  surface  waves  and  current*, 
we  substibute  Equations  (2.4.6)  into  (2.4.3)  obtaining, 

‘^t  7  +  gh  a  -  ^  +  gxJ 

at  z  =  x+h  .  (2.4.8) 

To  simplify  this  equation,  we  could  implement  (2.4.7a);  however. 
Equation  (2.4.7a)  is  applicable  only  at  z  =  x*  We  therefore 
expand  the  right  hand  side  of  this  equation  about  the  unperturbed 

*  We  will  use  the  term  current  to  describe  the  gradient  of  the 
long  wave  \’elocity  potential  regardless  of  the  generating 
mechanism. 
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surface.  We  obtain  the  following  expressions  at  the  unperturbed 
surface 


V'F  I 

' z=h+X 


V4'  +  h  1^  (VH')  +  .  .  .  j  I 


and 


4' 


t 


z=h+x 


+  ... 


z=X 


in  which  we  will  keep  only  terms  linear  in  the  surface  wave 
height  (h) .  Substituting  these  expansions  into  Equation  (2.4.8), 
we  obtain  for  the  interaction  equation 


+  2  V())*V4i  +  l^g  +  +  VV*V(()  + 

V4/  .yv 

|h  +  V4'*V(}) 

=  0 

which  to  terms  linear  in  (Ji  and  h  is 

(|)^  +  V4'-V<|)  T  j^g  +  V4'-VV  +  ||jh 

=  0 

at  z=x  . 

(2.4.9) 

i 

1 

1  1 

We  have  introduced  the  notation  in  the 

above 

equations 

that  1 

f.  ' 

the  surface  current  may  be  written  as 

) 

I 

1 

V4'  =  e^U(x,t)  +  e^V(x,t; 

! 

(2.4.10) 

where  U(x,t)  and  V(x,t)  are  the  horizontal  and  vertical  com¬ 
ponents  of  the  current,  respectively. 

1  "1,. 
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Equation  (2.4.9)  is  seen  to  retain  a  linear  inter¬ 
action  between  the  current  given  by  Equation  (2.4.10)  and  the 
surface  waves.  A  second  aspect  worthy  of  note  is  the  modifi¬ 
cation  in  the  gravitational  acceleration  due  to  the  finite 
value  of  the  vertical  current  component  which  induces  an  addi¬ 
tional  vertical  acceleration. 

We  may  write  Equation  (2.4.9)  in  the  alternate  form 

°  atz  =  x. 

(2.4.11) 

The  generating  mechanism  for  the  current  determines  the  relative 
magnitude  of  the  unperturbed  horizontal  and  vertical  velocity 
components.  For  example,  if  the  current  is  produced  by  an 
internal  wave,  then|u| »|v| so  that  only  the  horizontal  velocity 
component  need  be  retained  in  Equation  (2.4.11).  If,  however, 
the  current  is  produced  by  a  swell  passing  under  the  region  of 
interest,  the  vertical  component  of  the  velocity  may  be  com¬ 
parable  to  the  horizontal.  In  this  case,  both  components  must 
be  considered.  The  particular  mechanism  being  investigated 
will  therefore  determine  the  form  of  the  equation  used  in  the 
analysis. 

To  determine  the  linear  form  of  the  kinematic  boundary 
condition  [Equation  (2.4.4)],  w«  encounter  the  same  difficulty 
as  with  (2.4.3),  that  is,  going  from  the  free  surface  to  the 


unperturbed  surface.  To  obtain  our  expression,  vve  expand 
the  right-hand  side  of  Equation  (2.4.4)  about  z=X  as  follows: 


94' 

9z 


z  =h+x 


lil  +  h  ^  + 


2=X 


and  using  the  definition  of  the  current  components  in  Equation 
(2.4.10)  obtain , 


94* 

9z 


<94'  ^  .  9V/ 

z=h+x 

j9z  ^  9zj 

z=X 


to  terms  linear  in  h.  Of  course,  we  must  also  expand  the 
appropriate  terms  on  the  left  side  of  Equation  (2.4.4)  and 
use  the  expansion  for  V4'.  The  inclusion  of  these  expansions 
and  using  Equation  (2.4.8)  to  simplify,  we  obtain 


+  V  (t)*V  +  V4'*V  +  VV*V 


) 


h 


^  +  h  -  V(J:-VX  -  VV-VX 

o  Z  0  Z 


or  to  terms  linear  in  ({>  and  h 


V4'*V 


_  9  ({) 

9z 


V4cVX 


at  z  =  X.  (2.4.12) 


We  may  use  the  incompressibility  condition  (V»u  =  0)  to 
9  V  9U 

replace  and  write  (2.4.12)  as 


“  H  "  at  z-X.  (2.4.13) 

If  we  restrict  our  analysis  to  a  single  horizontal  dimension, 
we  may  write  (2.4.13)  as 


which  is  the  vertical  velocity  of  the  surface  waves.  We  note 
that  this  equation  was  assumed  to  have  the  form 


in  PD-7 2-023  at  the  perturbed  surface.  The  effect  of  calculating 
at  the  unperturbed  surface  is  to  modify  the  equation  by  the 
gradient  of  the  horizontal  surface  current  and  the  horizontal 
variation  in  the  unperturbed  surface  height. 

We  may  replace  the  horizontal  gradient  of  the  unperturbed 
surface  in  Equation  (2.4.14)  by  taking  the  horizontal  derivative 
of  Equation  (2.4.8)  to  obtain 


dX 

5T? 


Also,  using  the  fact  that  the  ocean  is  irrotational ,  i.e., 
V  X  u  >  0,  we  may  replace  with  so  that 
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r 


j 


3x 


1  ^ 
g  dt 


3x 


+  V 


(2.4.15) 


Substituting  Equation  (2.4.15) 


3t  3x  3z  3x  1 


into  Equation  (2.4.14) 

ii  +  i  ^ 

3z  g  3x  dt 


yields 


(2.4.16) 


as  our  second  linearized  dynamic  equation  for  the  surface  waves. 
Equation  (2.4.16)  has  the  same  structure  as  the  expression 
obtained  by  Zachariasen  (1972) ,  however,  it  is  evaluated  at 
z=x(x,t)  and  not  z=0  as  are  his  expressions. 

Milder  (1972)  has  shown  that  the  simplicity  of  these 
equations  can  be  maintained  by  noting  that  the  right-hand 
side  of  Eq.  (2.4.14)  can  be  written  as 


Ik  -  3X  -  9s  3<j) 
3Z  3x  3x  ~  3x  3n 


(2.4.17) 


where  s  is  the  distance  along  the  surface  z=x(x,t)  and  —  is 

dn 

the  derivative  normal  to  that  surface.  Since  the  velocity 
potential  satisfies  the  two  dimensional  Laplacian, 


Milder  also  shows  that 


3s  3(j)  _  _ .  3^ 
Jx”  5n  ”  ^  Sx” 


(2.4.18) 


We  premultiply  Eq.  (2.4.11)  by  the  operator 


^  +  M  i_  +  iu 

5T  3x  3x 
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I 


r 


and  use  Eqs.  (2.4.16)  and  (2.4.18)  to  obtain 
\3t  ”  3x  3xJ)3t  ^  3xj^''  3x 


(2.4.19) 


where  we  have  neglected  the  vertical  velocity  V  compared  with 
U  and  *  g  +  We  obtain  the  term  by  term  expansion  of 


Eq.  (2.4.19) 


^  +  u2  -ii  + 


A  H  "  ^^eff 


1  _L  +  A 

3x  9x  at 


(})  =  0 


(2.4.20) 


which  is  our  linearized  dyn2unic  equation  for  the  surface 
velocity  potential  at  the  unperturbed  surface  z=x(x,t). 


2.4B  Linear  Energy  Transfer  Equation 


To  investigate  Equation  (2.4.20),  we  assume  a  solution 


of  the  form 


A(x,t)  e 


ir(x,t) 


(2.4.21) 


where  the  frequency  and  wavenumbers  for  the  wavetrain  are  de¬ 


fined  by 


k  = 


and  (1) 


3r 

Jt  • 


(2.4.22) 


Substituting  Equation  (2.4.21)  into  Equation  (2.4.20)  and 
using  the  definition  of  the  frequency  and  wavenumber,  we  obtain 


9e£f  =  5  ^  *  |j}{^  +  « 


A  (2.4.23) 


when  |u|  >>  |v| ,  for  the  real  part  of  Equation  (2.4.20).  The 
imaginary  part  of  (2.4.20)  yields  the  equation 


J 


-  |tt  ^  u  geff 


0. 


(2.4.24 ) 


i£ 

3x 

All  the  terms  on  the  right-'iand  side  of  Equation  (2.4.23)  are 
of  second  order  and  in  the  WKB  approximation  we  may  write  the 
dispersj.on  relation 

(w  -  kU)^  -  k  «0  .  (2.4.25) 

If  we  introduce  the  following  definitions  of  a  velocity  (c) 
and  energy  (E)  , 

^  (2.4.26) 

and 

E  =  kA^  =  kUl^  (2.4.27)  I 

1 

I 

we  obtain  the  energy  conservation  equation 

It  ItT  [(“  +  f)'^]  '0  (2.4.28) 


from  (2.4.24)  using  (2.4.25).  Equation  (2.4.28)  may  bo 
written  in  the  alternate  form 


2.5  Wave  Packet-Surface  Current  Interaction 


In  this  section,  we  restrict  our  attention  to  !:he  case 
of  a  surface  current  whose  horizontal  velocity  component  is 
much  greater  than  the  vertical  component;  e.g.  ,  a  current 
generated  by  a  non-dispersive  internal  wave.  If  we  consider 
a  surface  wave  packet  moving  into  a  region  of  non-uniform  but 
temporally-steady  current,  we  can  write  the  average  value  of 
the  square  of  the  wave  packet  amplitude  in  terms  of  the  energy 
density  F(k)  aa 


a 


(2.5.1) 


The  wave  packet  is  considered  to  be  localized  in  a  region  Ax 
of  physical  space  and  must  therefore  be  concentrated  in  a  region 
Ak  of  wavenumber  space.  The  integral  Equation  (2.5.1)  is 
restricted  to  this  region  of  wavenumber  space  (Ak) .  If  the 
width  of  the  wave  packet  is  narrow;  i.e.,  Ak  <<  k,  then  the 
initial  wave  train  is  nearly  monochromatic  and  we  may  write 

~7 

a  ~  F(k)  Ak  .  (2.5.2) 


The  quantity  P(k)  in  Equation  (2.5.2)  is  the  average  energy 
density  of  ti^e  wave  packet  in  k-space.  We  will  drop  the  bar 
for  convenience  in  the  following  discussion. 


Let  us  consider  the  case  of  a  steady  horizontal  current 
U  =  U  (x) .  The  total  time  variation  (following  the  wave  packet) 


k 

I  » 
^  , 


Y 

\ 


.1 


in  the  squared  eunplitude  may  be  written  using  Equation  (2.5.2) 
as 


dt 


2 


a 


Ak||+F 


dAk 

dt 


(2.5.3) 


It  is  clear  from  Equation  (2.5.3)  that  as  the  wave  packet 
traverses  the  current  region  its  width  changes  markedly.  If 
we  define  the  packet  width  as  Ak  =  kj^-k2»  then 


dAk 

3t” 


dk^  dk^ 

3t“  ■  3t“ 


(2.5.4) 


but  we  know  from  the  dispersion  relation  of  the  preceding  sec¬ 
tion  that 


dk 

art  = 


-k 


3U 


so  that  the  time  variation  in  the  packet  width  given  by  Equa¬ 
tion  (2.5.4)  can  be  written  as 


dAk 

3^ 


-Ak 


3U  /4U  +  3c^ 
Uu  +  Ic) 


(2.5.5) 


Using  this  equation  in  Equation  (2.5.3)  and  dividing  the  re¬ 
sulting  expression  by  Equation  (2.5.2)  yields 


1  da^  «  i  ^  /4U  t  3c\ 

“  F  dt  "  Tx  V4U  +  2c  J 

cl 


(2.5.6) 


as  the  percent  variation  in  the  packet  amplitude  as  it  moves 
into  a  region  of  non-uniform  current. 
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2.5A  Wave  Packet  Distortion 


An  alternative  expression  to  Equation  (2.5.6)  may  be 
obtained  by  using  the  energy  conservation  equation 

It  H  ^  k  [(“  *  1)°=]  *  0 

of  the  preceding  section.  For  a  steady  non-uniform  current, 
we  have 

c  =  ^gA 

and 

''G  =  “^1 

for  the  phase  velocity  in  the  absence  of  the  current  and  group 
velocities  of  the  surface  waves,  respectively.  For  these  con¬ 
ditions  Longuet-Higgins  (1961,  1964)  among  others,  has  shown 
that  Equation  (2.5.7)  has  the  solution 

c v„E  =  const.  (2.5.8) 

\J 

To  express  Equation  (2.5.8)  in  terms  of  the  wavepacket 
amplitude,  we  introduce  the  quantities 

=  Ae^^  and  h  =  ae^^  (2.5.9) 

where  the  wavenumber  and  frequency  is  defined  by  Equation 
(2.4.23).  The  surface  elevation  and  velocity  potential  defined 


by  Equation  (2.5.9)  clearly  obey  Equation  (2.4.11), 


1_ 

3t 


+ 


U 


4)  +  gh  *  0 


(2.5.10) 


where  we  have  neglected  the  effects  of  the  finite  vertical 
velocity  V.  Using  Equation  (2.5.9)  in  (2.5.10),  we  obtain 


1 

X 


(w  -  ku)  +  g  I 


0 


(2.5.11) 


which  using  the  dispersion  relation  Equation  (2.4.26)  becomes 

"k  (fr  ”  lx)  A  i  -  I  .  (2.5.12) 


If  the  fractional  change  in  the  amplitude  is  much  smaller 
than  the  acceleration  of  gravity 


1  4A 

X  ar 


«  g 


in  the  WKB  limit,  we  obtain  from  Equation  (2.5.12) 


A  w  c  a 


(2.5.13) 


We  may  use  Equation  (2.5.13)  in  the  energy  conservation 
Equation  (2.5.8)  and  using  Equation  (2.4.28)  obtain 
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*\ 


,3  2  ^  ^ 

kc  Vq  a  =  constant 

g  so  that 


T 


c  v_  a  =  constant 


1 


: 


Returning  now  to  our  original  problem,  i.e 
a  second  expression  for  the  time  variation  in  the 
amplitude,  we  take  the  derivative  of  Equation  (2. 
follows 


d 

3t 


c 


0 


For  a  steady  current,  i.e.,  -g-  »  0,  and  using 

dc  _  _  c  dk  _  c  3U 
3t  ~  2Tc  3r  “  7  3x 


we  obtain 


1  da^ 


1  + 


1  U  t  c 
^  u  + 


3U 
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as  our  second  expression.  Equating  Equation  (2.5. 
Equation  (2.5.19)  yields 


(2.5.14) 

.  ,  finding 
wave  packet 
.14)  as 

(2.5.18) 


(2.5.19) 

6)  with 


I 

I 


1  dF 
F  3t  “ 


&U  /  4U  +  3c  \ 
3x  V  4U  +  2c  / 


1  U  +  c 


3x 


or 


1  dF 
F  at 


+  ,  1  3U 
^  7  97° 


0 


(2.5.20) 


which  accounts  for  the  non— dispersive  effects  associated  with 
compression  of  the  wave  packet  . 


2.53  The  Effect  of  Capillary  Waves  on  Wave  Packets 

The  effect  of  capillary  waves  may  be  included  in  the 
radiation  stress  obtained  in  the  preceding  section  by  modi¬ 
fying  the  appropriate  terms  in  Eq.  (2.5.7).  As  will  be 
discussed  in  the  following  section,  if  we  define  the  para¬ 
meter  Y*T/pg  where  T  is  the  surface  tension  and  p  the 
density  of  the  ocean,  the  frequency  of  surface  waves  becomes, 


0) 


gk (1  +  Yk  ) 


The  phase  velocity  c  is 


0) 


c  =  ^  =  /g/k  |/l  +  Yk^ 

and  group  velocity  becomes 


4  ^ 


1  +  3Yk' 


(2.5.21) 


(2.5.22) 


do) 

ai^ 


(2.5.23) 


which  becomes  c/2  as 

We  can  use  these  definitions  in  the  conservation 
expression  Eq.  (2.5.18)  to  yield 


d 

at 


{ 


c(U  +  c  ) 
g 


0 


(2.5.24) 


For  a  steady  current,  i.e., 


3U 

3t 


0  ,  and  using 


(2.5.25) 


(2.5.26) 


where 


dt 


(U+Cg 


we  obtain 


1  ^  _ 
^  dt  ~ 


from  Eq.  (2.5.24) 
1  +  1  |(2+U/Cg)  - 


12Yk^  ^ 
l+3Yk‘ 


(2.5.27) 


in  place  of  Eq.  (2.5.19).  Substituting  Eq.  (2.5.27)  into 
Eq.  (2.5.6)  we  arrive  at 


(2.5.28) 


1  dF  .  1 
F  dt  I 


1  -  Yk^ 

1  Ty? 


3x 


0 


as  the  modified  radiation  stress  term  due  to  the  action  of 
capillary  waves. 


2. 5C  Non-Linear  Amplitude  Effects  on  Effective  Group  Velocity 

The  dispersion  relation  used  in  the  last  section  was 
of  course  that  given  for  a  linear  gravity  wave  modified  to 
include  capillary  waves.  We  now  wish  to  model  the  change 
in  the  dispersion  relation  due  to  finite  amplitude  effects 
by  incorporating  in  our  equations  a  modified  group  velocity. 

A  more  systematic  approach  to  this  problem  would  be  to  modify 
the  energy  balance  equation  introduced  in  Section  2.4  to 
include  the  finite  amplitude  effects.  This  approach  is  presently 
being  explored  and  will  be  reported  at  a  later  time. 

To  second  order  the  dispersion  relation  of  ar  isolated 
non-linear  gravity  wave,  i.e.,  Stokes  wave,  is 

CO  -  kc  »  */glc  (1  +  k^a^)^  (2.5.29) 


where  ka  is  the  slope  of  the  wave.  The  group  velocity  of 
this  wave  is 

Cg  -  -  c^°’  (l  +  2  kV)  / yJTTk^  (2.5.30) 
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wher''  is  the  linear  group  velocity  ^s/g/k  or  that  given 

g 

by  Eq.  (2.5.23).  To  generalize  Eq.  (2.5.30)  to  the  situation 
where  there  is  a  spectrum  of  waves,  we  introduce  the  notion 
of  a  correlation  interval.  We  assume  that  waves  in  the 
interval 


k-^<k<k  + 


Ak 

T 


are  correlated,  "aves  outside  this  interval  are  assumed 
uncorrelated  and  not  to  alter  the  wave  phases  in  any  syste¬ 
matic  way.  Thus  we  write  the  dispersion  relation 
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If  the  spectrum  is  sharply  peaked  we  write 


Cg  =  ^  I  •^^F(k)Ak} 


(2.5.33) 


The  model  we  use  in  our  calculation  parameterizes  the  effect 
of  Ak  in  Eq.  (2.5.33)  so  that  we  have 


=9  =  =9°’  (l  +  >  ■>'^F(k)) 


(2.5.34) 


where  A  is  a  parameter  on  the  order  of  unity. 
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A _ -  ^ 


A 


>■ 


2-6  Model  for  Generation  of  Waves  by  Wind;  Effects  Independent 
of  Wave  Slope 


In  this  section  we  wish  to  outline  the  initial  growth 
of  surface  waves  due  to  the  statistical  fluctuation  in  air 
pressure  due  to  turbulence,  eddies,  wind  gusts,  etc.  To  con¬ 
struct  this  model,  we  consider  the  linear  form  of  the  inter¬ 
action  equation  for  infinitesmal  surface  gravity  waves  dis¬ 
cussed  in  Section  2.4.  We  modify  the  previous  equation  by 
the  inclusion  of  the  incremental  pressure  with  respect  to 
ambient  (Ap) ,  that  is 


;^+gh+-Ap=0  atz=0  (2.6, 

and 

dh  3  4 

at  “  at  z  =  0  .  (2.6. 

We  note  that  the  above  equations  are  taken  about  the  z  =  0 
position,  which  implies  that  there  is  no  ambient  current 
present,  i.e.,  u  =  V  =  0  from  Section  2.4. 

The  expression  for  the  pressure  increment  may  be 
obtained  by  considering  the  interface  between  two  fluids. 
Consider  the  vertical  forces  acting  on  a  strip  of  surface 
of  width  6x.  If  p  is  the  pressure  just  below  the  water 
surface,  p  the  pressure  of  the  air  and  T  the  surface  tension 


then 


(P  -  P3)6x  V  s[t 

is  the  force  balance  acting  on  the  strip  6x.  The  pressure 
is  therefore  given  by 


The  air  pressure  p  can  be  modeled  by  a  sum  of  these  terms; 

A 

(i)  the  direct  coupling  of  the  wind  to  the  water,  (ii)  the 
local  statistical  variation  in  the  surface  pressure,  and  (iii) 
the  surrounding  or  ambient  pressure  at  the  surface. 


(2.6.3) 


P.  -  <>'‘11*''  +  Panbient  '  <2.6.5) 

The  first  term  on  the  right  of  Equation  (2.6.5)  gives  the  in- 
phase  pressure  variation  at  the  surface  due  to  the  action  of 
wind.  This  term  will  be  discussed  more  fully  in  Section  2.7. 
Substituting  Equation  (2.6.5)  into  Equation  (2.6.1)  yields 


d(()  .  -v  .  3h  T 

3^  +  gh  +  a  -  - 


T 


at  z  =  0 


(2.6.6) 


To  obtain  a  solution  to  Equation  (2.6.6)  and  (2.6.1), 
we  decompose  the  vector  potential  and  surface  waves  into  their 
fourier  components  and  thereby  obtain  a  modal  description  of 
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our  linear  problem.  The  surface  elevation  is 


h 


ikx 
e _ 

/L 


(2.6.7) 


and  since  (j)  obeys  Laplace's  equation  [V  (j)  =  0]  the  vector 
potential  is 


<P 


-  E  K  (t) 


ikx 
e _ 

/L 


(2.6.8) 


Usinq  Equation  (2.6.2)  we  obtain  the  relation  between  the 
mode  amplitudes 


and  because  our  dynamic  equation  is  taken  at  z  =  0,  we  obtain 

hk  +  jgJ'  +  \  =  "I  (2.6.10) 

when  Equations  (2.6.7)  -  (2.6.9)  are  substituted  into  Equation 
(2.6.6).  The  quantity  P(k,t)  in  Equation  (2.6.10)  is  the 
statistical  fluctuation  of  the  pressure  in  k-space  and  y  =  T/p. 

The  structure  of  Equation  (2.6.10)  is  clearly  that  of 
a  driven  harmonic  oscillator 

^k  “k  \  =  ~  7  .  (2.6.11) 


The  driving  force  in  this  case  is  the  statistical  'unction 
P(k,t).  This  function  is  uncorrelated  with  the  gravity  wave 
amplitudes  [hj^  (t)]  and  acts  to  randomize  the  growth  of  the 
mode  amplitudes. 

The  solution  to  Equation  (2.6.11)  in  the  homogeneous 
case,  with  an  initial  amplitude  and  slope  of  zero,  i.e., 
hj^(t  =  0)  =  hj^(t  =0)  =  0,  is 


(2.6.12) 


The  solution  to  the  inhomogeneous  equation  can,  of  course,  be 
written  in  terms  of  Equation  (2.6.12)  as 


k 


P(k,t)  sin  )  dt' 


(2.6.13) 


which  determines  the  statistical  character  of  the  mode  amplitudes. 
We  can  smooth  over  these  fluctuations  by  means  of  an  ensemble 
average  and  concern  ourselves  with  the  spectrum  of  surface  dis¬ 
placements.  Using  Equation  (2.5.2),  we  obtain 


F(k,t) 


<hj^(t)  h*j^(t)> 

Alt 


(2.6.14) 


as  the  quantity  of  interest,  rather  than  the  mode  amplitudes 
themselves.  The  brackets  in  Equation  (2.6.14)  denote  an 
ensemble  average  and  Ak  is  the  element  of  wavenumber  space 


occupied  by  the  wave  packet  of  interest. 

Substituting  Equation  (2.6.13)  into  (2.6.14),  we 
obtain  for  the  surface  spectrum 

X  sin  {i(t-t')  sin  f2*(t-t")  (2.6.15) 

where  <P  (k  ,t ' ) P  (k , t" )>  is  the  spectrum  of  the  pressure  fluctu¬ 
ations.  If  we  approximate  the  comp] ex  frequency  as 

~  \/a)^  +  ik^a  =  ia/c^  »  tOj^d  +  ia/2c^  ) 

(2.6.16) 

and  assume  that  the  statistical  pressure  field  is  stationary, 
Equation  (2.6.15)  can  be  simplified.  It  is  not  necessary  to 
reproduce  here  the  detailed  reduction  of  this  equation  [see 
Phillips  (1966) ,  p.  85] ;  we  merely  note  that  when  the  time 
interval  is  greater  than  the  time  over  which  pressure  components 
of  wavenumber  k  are  correlated.  Equation  (2.6.15)  can  be  written 
in  the  form 

F(k,t)  = 

p  c^ 
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The  quantity  M  =  a/c  is  the  coupling  parameter  between  the 
wind  and  water. 

In  the  regime  where  the  hyperbolic  sine  may  be  expanded, 
that  is,  for  time  t  <<  ^  ,  Equation  (2.6.17)  may  be  written 
as 

F(k,t)  =  Q(k,a))  t  .  (2.6.18) 

pc 

Alternatively,  we  may  write  for  the  rate  of  growth  of  the 
energy  spectrum  close  to  the  time  t  =  0,  i.e.,  the  initial  wave 
growth 


3F 

Tt 


Q(k,a3) 

P  c 


(2.6.19) 


The  functional  form  we  use  for  Q(k,u))  is  that  proposed  by 
Barnett  (1968)  which  is  based  on  the  work  of  Priestly  (1965) 


Q(k,a))  = 


6.13xl0”V 

2  2 
TT  U 


V, 


k 


sin  0 


+  (k  cos  0  -  H) 


(2.6.20) 


where  the  following  definitions  are  used; 
p  is  the  density  of  water  (kg/m^) 


c  is  the  surface  wave  phase  velocity  (=  w/k) 


t 

t 

I 

I 

9  is  the  angle  between  the  wave  and  wind  vectors 
W  is  the  wind  speed  (m/sec) 

H  is  u)/W  (1/m) 
is  0.33 

V2  i3  0.52 

These  values  of  and  V2  are  taken  from  Collins  (1972)  . 


2 . 7  Model  for  the  Generation  of  Waves  by  Wind;  Wave  Slope 
Dependent  Effects 

The  wind  blowing  over  the  surface  of  the  ocean  is  the 
primary  mechanism  for  the  generation  of  waves.  The  wind  sets 
up  a  pressure  distribution  at  the  ocean  surface  and  it  is  this 
variation  in  pressure  which  stimulates  the  growth  of  waves . 
Phillips  divides  these  pressure  fluctuations  into  two  broad 
categories  :  those  induced  by  the  coupling  of  the  wind  to  the 
irregular  ocean  surface  and  those  produced  by  turbulent  eddies 
in  the  air  flow.  The  first  effect  is  treated  as  a  direct 
coupling  of  the  ocean  wave  field  to  the  average  wind  field. 

The  second  effect  is  treated  as  a  random  disturbance  on  the 
ocean  wave  amplitudes  of  different  wavenumbers.  Both  these 
terms  feed  energy  into  the  surface  waves;  the  direct  term 
into  a  narrow  band  around  k,  which  corresponds  to  the  matching 
of  the  average  wind  speed  and  phase  velocity  of  the  surface 
wave,  and  the  turbulent  term  into  a  broad  spectral  band  around 
JT. 

Assume  that  the  strength  of  the  coupling  between  the 
wind  and  water  is  given  by  the  complex  parameters  (v  +  iy) •  It 
has  been  determined  that  the  real  part  of  the  parameter  (v)  is 
generally  negative  and  models  the  dissipation  of  energy  due  to 
the  formation  of  spray  at  wave  crests.  The  negative  value  of 
V  produces  a  decrease  in  pressure  (suction)  at  a  wave  crest  and 
an  increase  in  pressure  at  the  trough,  thereby  aiding  in  the 
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deformation  of  the  wave.  This  model  of  spindrift  has  been 
included  in  principle  within  the  wave-breaking  terms  and  will 
not  be  further  considered  here.  The  imaginary  part  of  the 
coupling  pareuneter  (Y)  models  the  rate  at  which  new  waves  are 
generated  in  the  ocean  wave  field  and  is  discussed  more  fully 
below. 

The  coupling  of  the  wind  to  water  waves  determined  by 
the  above  pareimeter  produces  a  shearing  force  at  the  ocean 
surface.  Phillips  (1966,  Section  4.j)  discusses  two  contri¬ 
butions  to  this  stress:  one,  a  resonant  contribution  from 
the  region  in  the  boundary  which  moves  at  the  same  speed 
as  the  wave  under  consideration  (the  matched  layer) ,  and 
the  other,  the  contribution  of  the  boundary  layer.  The 
matched  layer  contributions  dominate  in  the  region  of  strong 
interaction  and  is  the  only  effect  we  consider  here.  For 
this  contribution,  Phillips  adapted  a  model  developed  by 
Miles  (1957)  for  quasi-laminar  flow.  In  this  model  the 
Reynolds  stress  is  written  in  the  form 


2,  7  2 
2YPwC  k 


9 


(2.7.1) 


t 


where  is  the  density  of  water  and  Y  is  the  imaginary  part 
of  the  coupling  parameter.  We  may  relate  the  mean  momentum 
flux  'l>(k)  to  the  mean  energy  flux  H'(k)  by  the  equation 


<I)(k)  =  — IS—  4'(k) 
a)(k) 


(2.7.2) 


so  that 


(if)., 


=  Y  F(k) 


wind 


(2.7.3) 


where  we  have  used  the  relation 


between  the  energy  density  and  energy  flux. 

The  coupling  parameter  y  is  modeled  by  considering 
the  mean  flow  over  a  (moving)  wavy  surface.  The  primary 
contribution  to  Y  is  from  the  region  where  the  wind  speed  is 
just  equal  to  the  wave  speed,  i.e.,  at  a  height  Phillips 

calculates  the  contribution  to  Y  from  this  matched  layer  to  be 


(2.7.5) 


« 


/ 
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The  equilibrium  wind  profile  is  approx iK=itely  of  the  form 


U(z)  -  c  *  ^  log  (2.7.6) 

o  V  m/ 

where  U*  =  is  the  wind  friction  velocity  and 

0.42  is  Karman's  constant.  With  this  substitution  into 
the  above  equation,  we  have  as  the  coupling  parameter 


(2.7.7) 


where 


B  =  (kz^) 


r 


log  u  exp(-kZj^u)  du 


(2.7.8) 


If  from  the  equilibrium  wind  profile  the  point  where 
U(z^)  =  0  is  taken  as  the  origin,  then  for  an  aerodynamically- 
rough  flow 


*  — ^2-  il^j  (2.7.9) 

A  2  ^2 

where  z^  %  B^U  /g  and  sr  1.1  x  10  as  suggested  by 
Charnock  (1955)  .  This  empirical  relation  sesiTiC  to  be  a  good 
representation  at  all  but  the  lowest  wind  speeds,  where  the 
flow  is  aerodynamically  smooth. 


2 

Miles  found  A  to  be  equal  to  n  and  F  ~  1/3  for  z  >  z  , 
m  ^  '  m 

and  that  these  values  are  insensitive  to  fluid  viscosity  in 
quasi-leiminar  flow.  In  the  present  calculation  then,  these 
terms  are  constants  in  the  above  equations.  In  Figure  4  the 
quantity  3  is  plotted  as  a  function  of  the  ratio  C/U*.  It 
is  clear  from  this  figure  that  the  coupling  of  the  wind  to  the 
waves  is  restricted  to  values  of  C/U^  S  20.  Waves  having 
phase  speeds  greater  than  20  U*  effectively  "out  run"  the 
wind  and  are  not  amplified  significantly. 

Measured  in  wave  periods,  the  growth  time  [T  =  (ajY)~^sec] 
is  given  by 


»  26.8  ;c/U*)^/3 
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FIGURE  4. 


Coupling  of  the  Wind  to  the  Ocean  Surface 
as  a  Function  of  C/U* 
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2 . 8  Wave  Cresting  Effects 

As  the  spectrvun  develops  into  a  fully  saturated  state, 
interactions  between  waves  having  widely  differing  wave-lengths 
alter  the  wave-breaking  phenomena  discussed  in  Section  2.2. 

Phillips  (1963)  has  identified  this  mechanism  as  being  primarily 
responsible  for  the  attenuation  of  long  gravity  waves  passing 
through  a  saturated  wave  region.  Here  as  the  long  wave-length 
waves  or  swells  pass  under  the  short  gravities,  the  wave  breaking 
at  the  crests  of  the  long  waves  is  enhanced  whereas  in  the  troughs 
the  waves  are  calmed.  Phillips  estimates  the  attenuation  rate 
of  the  long  waves  by  considering  the  time  dependent  exchange  of 
energy  between  the  short  and  long  waves .  In  the  absence  of 
breaking  or  wind  input,  this  exchange,  which  results  from  the 
radiation  stress  term,  is  oscillatory  and  averages  to  zero.  In 
the  presence  of  wave  breaking,  there  is  a  net  transfer  of  energy 
from  the  long  waves  to  the  short  waves,  which  in  a  steady  sea 
is  balanced  by  the  wind  input. 

The  seune  phenomena  should  be  responsible  for  a  calming 
effect  of  long  waves  on  short  waves.  Here,  when  the  sea  has  de¬ 
veloped  to  the  point  where  long  waves  have  appreciable  amplitude, 
there  is  a  net  enhancement  of  breaking  of  the  short  gravities 
over  that  which  occurred  before  the  long  waves  had  developed. 

The  net  result  then  is  that  the  mean  wave  height  or  slope  of 
the  short  waves  in  the  saturated  region  tends  to  decrease  as 
the  long  wave-length  portion  of  the  spectrum  develops. 


In  this  section  we  construct  models  for  both  the.sc. 


effects . 


2.8A  Ef ':ect  of  Long  Waves  on  Short  Waves 

We  consider  the  alteration  of  the  wave  breaking  term 
induced  by  a  long  wave-length  wave.  The  wave-breaking  term 
has  the  form  (from  Section  2.2) 


breaking 


(-0^  /  k^F) 


(2.8.1) 


We  may  write  F  in  the  form 

F(k,x)  =  f(k)  +  Fj^(k,x)  (2.8.2) 

where  F(k)  is  the  mean  value  of  the  wave  height  spectrum 
averaged  over  the  long  waves  and  Fj^(k,x)  is  the  perturbation 
of  the  short  wave-length  spectrum  induced  by  the  long  waves. 

To  estimate  the  effect  of  the  long  waves,  we  expand  the 
right-hand  side  of  Equation  (2.8.1)  to  second  order  in  F^^.  It  is 
convenient  for  this  purpose  to  define  a  function  G(y)  according  to 


G(Y) 


y3/2  ^-1/Y 


(2.8.3) 


Then 


(—) 

V^^/breaking 


Qq 

k^/? 


G(kV0o) 


(2.8.4) 


I 


and  approximately 


© 


breaking 


k^/iT 


where 


= 


(2.8.5) 


(2.8.6) 


=  k^F/02 


(2.8.7) 


YGy/G 


(1  +  I  y)/Y 


(2.8.8) 


(1  +  Y  +  I  Y^)/Y^ 


(2.8.9) 


For  linear  gravity  waves  riding  on  a  steady  current  U (x) ,  we 
indicated  in  Sections  (2.4)  and  (2.5)  that  the  following  con¬ 
servation  rules  apply. 


k  ( /g/k  +  U)  =  constant 


a^  /g/k 


(u  = 


constant 


(2.8.10) 


(2.8.11) 


Thus  a  small  change  in  current  6U  produces  an  amplitude 
change  6a  according  to 


+  3U 


(/gTk +2U) 


a 


r 


(2.8.12) 


■7^ 


I 


which  can  also  be  obtained  from  Equation  (2.5.19)  by  equating 
2  .  /  c  \  3a^ 

6a  with  (U  +  j  j  •  6t.  The  mean  current  U  in  this  case  is  to 
be  taken  as  the  phase  velocity  of  the  swell  or  long  wavelength 

i 

gravity  U  ~ /g/k  .  Since  most  of  the  long  wavelength  waves 

s 

are  moving  much  faster  than  the  short  gravities,  we  can  take 
U  >>  /g/k  and  write 


6a  3  6U 

—  -  ~T  — 


(2.8.13) 


Now 


using  the  relation  ~  2a6a  and  F  ~  a  ,  we  find 


3 

I 


u 


(2.8.14) 


For  a  long  wavelength  gravity  wave,  U  is  the  phase  velocity 


(  /g/k  )  and  6U  is  the  orbital  velocity 
s 


6u  =  (jj  a„  cos  (k„x  -  w„t) 
s  s  s  s 


(2.8.15) 


Thus 


Y 


TT  a  k  cos  (k  X  -  w  t) 
2  s  s  s  s 


(2.8.16) 


Averaged  over  the  period  of  the  long  wave ,  the  self-breaking 
term  becomes 


(If) 


breaking  k'^/if  \0 


1  4. 

l+Y+|-Y^ 

9  2.2 

X  T 

y2 

8  s  s 

(2.8.17) 
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This  calculation  assumed  a  single  long  wavelength  wave  of 

amplitude  a  .  For  a  spectrum  of  waves ,  it  is  appropriate  to 
s 

2  2 

replace  the  term  a  k  by 

S  8 


2  2 
ajk*^ 
s  s 


/ 


F(k)  k  dk 


(2.8.18) 


since  only  waves  having  wavelengths  greater  than  that  of  the 
gravity  wave  will  contribute  to  the  wave-cresting  effect. 

The  final  form  of  the  wave-breaking  expression  now 


becomes 


,  -eu0 
3F(k)  _  '^o 


at 


k^/TT 


G(Y) 


1  +  w 

3  2~ 
i+Y+Jy^ 

k  1 

1  F(k)k^dk 

^  ^ 

Y^ 

1 

with 


(2.8.19) 


3  2 

y  =  k  f/g;: 

o 


The  term  in  the  curly  brackets  is  the  enhancement  of  the 
breaking  of  short  waves  due  to  the  presence  of  the  long 
waves . 

When  the  spectrum  is  saturated  above  a  wavenumber  k  . 

^  min 

(i.e.,  Fk^  ~  constant  for  k  >  k  .  ) ,  we  may  approximate 

min 

t 

F(k)  k^dk --  Fk3  An  k/k 

mm 
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Thus  the  wave-breaking  term  takes  the  approximate  form  for 
a  partially  saturated  sea 


1 

^breaking 


F  9t 


{-OqA^F) 


Since  the  magnitude  of  the  enhancement  term  is  strongly 
dependent  on  the  choice  of  the  parameter  0^,  careful  comparison 
with  data  will  be  required. to  arrive  at  a  useful  model. 


2 . 8B  Effect  of  Short  Waves  on  Long  Waves 

As  pointed  out  in  the  introduction  to  this  section, 
Phillips  has  constructed  a  model  of  the  interaction  between 
swells  and  short  gravity  waves  in  which  the  short  waves  act 
as  a  damping  mechanism  for  the  longer  swells.  The  notion 
employed  is  that  a  swell  passing  through  a  region  of  saturated 
sea  will  transfer  energy  to  the  gravity  waves  by  means  of  the 
"radiation  stress."  This  transfer  enhances  the  breaking  of 
surface  waves  at  the  crests  of  the  swell  and  calms  the  surface 
in  the  swells'  trough.  During  this  process,  the  wind  con¬ 
tinually  feeds  energy  into  the  gravity  waves,  so  that  although 
energy  is  lost  at  the  swell  crest,  the  wind  builds  the  waves 
back  up  after  the  breaking. 

Following  Phillips  (1966,  Section  4.9)  we  use  the 
energy  balance  equation  from  Section  2.4: 


T 


3E  ^  3 

■5t  ^ 


[(U  +  |)E] 


Q((jl)E 


(2.8.20) 


where  U  is  the  surface  current  produced  by  the  swells  and  the 
right-hand  side  of  Equation  (2.8.20)  is  the  energy  fed  into 
the  surface  waves  by  the  wind.  We  recall  that  the  phase  velocity 
in  Equation  (2.8.20)  is  given  by 


=  -  ^(9  t  a^)A 


(2.8.21) 


since  the  vertical  acceleration  of  the  swell  is  finite.  Because 

the  surface  pattern  is  progressing  with  the  phase  velocity  of 

the  swell  (c^)  we  may  replace  with  -  ^  in  Equation  (2.8.20). 

s 

To  determine  the  effect  of  the  short  waves  on  the  swell,  we 
need  only  consider  t.ie  interaction  between  the  surface  waves 
and  current. 

During  a  single  cycle  of  the  swell,  the  variation  in  the 
energy  density  E  of  the  short  waves  is  small  relative  to  the 
total  energy  density.  We  can  therefore  approximate  the  inter¬ 
action  term  from  Equation  (2.8.20)  by 


3 

Tt 


3  UE 
I  ~ 


s 


(2.8.22) 


which  is  correct  to  second  order.  We  can  integrate  Equation 
(2.8.22)  to  find  the  net  energy  transfer  from  the  swell  to  the 
short  waves  over  a  cycle,  that  is 
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=  If 


(2.8.23: 


where  6E„  is  the  change  in  energy  of  the  short  waves  due  to 
the  direct  action  of  the  wind. 

If  T  is  the  tiTne  interval  necessary  for  successive 
crests  of  the  swell  to  overtake  a  packet  of  surface  gravity 
waves,  we  may  write  for  the  change  in  energy  of  the  swell 
during  this  time  to  be 


dE(kg) 


(2.8.24) 


During  this  same  time  interval,  we  can  write  for  the  energy 
fed  into  the  surface  waves 


6E^  =  awET 


(2.8.25) 


where  E  is  the  average  energy  of  the  surface  waves.  Using 
Equations  (2.8.25),  (2.8.24)  and  (2.8.23),  we  can  write  for 
the  change  in  swell  energy 


dE  (ks)  _  3  U  .  X  .  r ..  . 

— —  2  c~ 


(2.8.26) 


over  a  cycle. 


We  write  the  velocity  potential  for  the  swell  in  one 


dimension  as 
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(2.8.27) 


X  =  -c^  cos  [k  X  -  oj  t] 
s  s  s  s 

so  that  the  current  is 

U  =  k  c  a  sin  [k  x  -  oj  t]  .  (2.8.28) 

s  s  s  s  s 


We  can  write  for  the  maximum  change  in  the  swell  energy,  using 
Equation  (2.8.28) 


dE(kg) 

“31“ 


-  »  k  a  a  (k)  w  (k)  E(k) 
^  s  s 


(2.8.29) 


which  is  the  result  derived  by  Phillips  (pg.  152) . 

If  we  now  wish  to  generalize  Equation  (2.8.23),  we 
could  consider  a  spectrum  of  surface  gravity  waves,  that  is, 
since  from  Equation  (2.5.2) 


E(k)  S3  pgF(k)  Ak 


for  a  nearly  monochromatic  wave  train  (Ak  <<  k) ,  we  might  make 
the  replacement 


a(k)  u)(k)  F(k)  Ak 


a(k)  a)(k)  F  (k)  dk 


(2.8.30) 


so  that  Equation  (2.8.29)  becomes 


dE(kg) 

dt 


-  I  X 


a(k)  w(k)  F(k)  dk 


(2.8.31) 


t 
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For  a  single  swell,  we  may  write  E(k  )  =  pg  a  for  the 

S  s 

energy  density  and  E (k  )  »  pgF  (k  )  so  that  Equation  (2.8.31) 

s  s 

becomes 


dF(kg) 

3t 


yiFikg) 


a  (k)  0)  (k)  F  (k)  dk 


(2.8.32) 


as  the  rate  of  dissipation  of  the  swell  energy  density. 


2.9  Wind  and  Wave  Induced  Surface  Currents 


In  a  fully  developed  wind  driven  sea,  the  turbulent 
stress  acting  on  the  ocean  surface  continuously  transfers 
momentum  from  the  atmosphere  into  the  water.  In  a  developing 
sea,  part  of  this  momentum  appears  as  increasing  wave  momen¬ 
tum  primarily  among  the  longer  developing  waves.  In  a  fully 
developed  sea,  however,  this  momentum  eventually  appears 
as  a  wind  driven  mean  current.  At  steady  state,  a  tur¬ 
bulent  boundary  layer  develops  in  the  water  near  the  sea 
surface.  The  turbulent  stress  created  by  the  atmosphere  is 


"a  =  Pa  ^*a 


(2.9.1) 


and  is  resisted  by  the  turbulent  stress  transferred  through 
the  water 


The  U*  is  the  friction  velocity;  the  remaining  notation  is 
obvious.  In  steady  state  x.  ~  t 

OL  W 

Using  a  similarity  argument  for  a  fluid  flowing 
over  a  smooth  surface,  t’^e  velocity  gradient  can  be  written 


where  k'^  ::  0.42  is  Karman's  constant.  We  may  therefore 
write  the  velocity  profiles  (see  Figiire  5) 


u  (z) 

a 


+  Uq  in  air 


(2.9.3) 


u  (z)  - - £n 

w 

o 


z  +z 


+  in  water 


(2.9.4) 


where  z^  is  a  roughness  length  and  Uj^  is  the  velocity  of 
the  water  surface.  We  may  relate  the  friction  velocities 
in  air  and  water  by  taking  the  steady  state  condition 
ss  ,  so  that 


/Pa\^  U. 

u  =  (  — )  u*  W  — - 
*w  \p../  *a  30 


(2.9.5) 


corresponding  to  typical  friction  velocities  in  water  of 
the  order  of  centimeters  per  second. 

We  take  the  total  current  to  be  zero  at  some  large 
water  depth  D.  The  surface  drift  velocity  is  then  given 
in  terms  of  D: 


"^o  vw/ 


D+z  '  ' 
w 

w 


(2.9.6) 
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FIGURE  5 


The  wind  and  Water  Shear  Profiles  Matched 
at  the  Oceaui  Surface  by  Equating  the  Shearing 
Stresses 


by  setting  u^(z)  =  0  in  Eq.  (2.9.4).  At  the  depth  |z| 
the  induced  current  is  given  by  Eq.  (2.9.4)  which,  using 
Eq.  (2.9.6),  becomes 


u  (z) 
w 


(2.9.7) 


A  proper  analysis  would  treat  the  wave  motion  in  a  shearing 
layer  in  detail  [such  an  approach  is  discussed  by  PhiDips 
(1973)].  Consistent  with  the  semi-empirical  approach  of  the 
current  model  and  in  view  of  the  very  weak  sensitivity  of 
the  result  to  our  assumptions  about  the  profile  we  adopt  a 
very  simple  approach.  We  assume  that  the  wave  motion  in  the 
presence  of  a  shearing  current  is  essentially  the  same  as 
that  of  a  wave  in  a  uniform  current.  The  value  of  this  ef¬ 
fective  uniform  current  should  depend  on  the  depth  to  which 
the  wave  orbital  velocities  penetrate.  To  find  the  effective 
mean  current  (u(k))  for  a  given  wave  (k)  ,  we  weight  this  pro¬ 
file  [Equation  (2.9.7)  with  a  suitable  power  of  the  amplitude 
of  the  wave  orbital  velocity. 
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where  oikD  >>  1  and  y  =  akz.  In  general,  the  quantity 

kz  is  extremely  small  and  the  lower  limit  of  the  integral 
may  be  replaced  by  zero,  yielding 


u„(k) 


•^o  WJ  V 


AnakD  +  y 


(2.9.9) 


where  y  is  Euler's  constant  {y  =  .5772...).  Substituting 
the  appropriate  values  for  the  parameters  in  Eq.  (2.9.9) 
yields 

u  (k)  »  0.0744  iln(1.73akD)  .  (2.9.10) 

Equation  (2.9.10)  gives  us  an  expression  for  the  average  drift 
current  felt  by  a  surface  wave  of  wavenumber  k.  This  weak  de¬ 
pendence  on  the  power  a  gives  us  confidence  in  this  method  of 
replacing  the  actual  sheared  boundary  layer  by  a  single  ef¬ 
fective  mean  current. 


07  07 
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2.10  Assembly  of  the  Transfer  Equation 


The  various  terms  modeled  in  the  preceding  sections 
may  now  be  assembled  into  the  overall  transfer  equation. 
Thus  we  have 


=  -  (c 

+  k 

+  ire 


g  +  UM 

3U'  ^ 

3x  3k 


^  3 

31? 


II 

3x 


(k^F)^ 
.  v/iT 


3 

3k 


(k^F) 


'] 


convection 

refraction 

wave-wave  interaction 


1  /l  -  Yk^ 
■  \l  +  Yk^ 


radiation  stress  with 
capillaries  plus  k- 
compression 


wave  breaking 


+  C2/glt 


6F 


matched  layer  wind 
generation 


-  4v^  k^F 


turbulent  and  viscous 
damping 


+ 


1.18  X 10 


-11 


7^ 


1  + 


1  -  c/w 


:  turbulent  wind 


(in  m  /sec) 


(2.10.1) 
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where  we  have  employed  the  following  definitions 


U'  =  U  +  u(k) 


from  Section  2.9 


(Y  =  k^F/0^) 


from  Section  2.5B,C 


from  Section  2.8 


H  =  in  meters"^  from  Section  2.6 

and 

i 

W  =  22.5  U*  in  meters/sec. 


The  transfer  equation  may  be  put  in  a  more  convenient  form 
by  introducing  a  number  of  dimensionless  variables: 

2 

wavenumber:  <  =  kU^/g 

distance;  C  =  x  g/U* 

time;  i  =  t  g/U* 

2 

surface  tension;  Y  =  ^  Y 

u; 

current;  u  =  U '/U* 

Y  =  F(k)kV0Q 


power  spectrum: 
angular  frequency; 


n  =  0)  U*/g 


In  this  system  of  units  we  have  the  relations 


and 


,  U* 

U'  +  c  =  - 

5  2/^ 


A  (  k)  +  2u/icj 


where 


A(k)  =  (1  +  c,y) 

a;?7  ' 

Using  these  relations  we  may  write  the  dimensionless 
form  of  the  transfer  equation  as. 


.u/7)||-2.3/ 


2  du  ^ 

3X  9< 


-  ono  7/2  3 


-  /IT  Y 


K 

3k 


-  —  ckY^'^^  (l+A)  exp(-l/Y) 
/tT 

+  C2K^6 ( K ) Y 


-  {8gv^/U*  )  K^/^  Y 


+  c. 


U*  \5.46  ^2.385 


0, 


1  +  ^1-1/22. 5/i^J 


-1 


(2.10.2) 


where 


s,  = 


.625  /^refV* 

c  ^.1.28  \  / 
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.0.28 


(22.5*  k) 
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and  C2=0.0119,  C2=1.89  and  is  variable. 

An  examination  of  Eq . (2 . lo . 2) shows  that,  in  the  absence 
of  viscous  and  surface  tension  effects,  the  value  of  the  wind 
speed  can  be  scaled  out  of  the  problem  by  expressing  the  data 
in  terms  of  the  reduced  variables  just  given.  The  general 
scaling  laws  are 


in  the  open  ocean,  and 


at  steady  state  in  coastal  waters  (x  is  the  distance  from  the 
lee  shore:  the  fetch) .  A  "fully  developed"  sei  develops  after 
a  long  time  in  the  open  ocean 

"k.  ’  (^  ' " ) 

To  show  self  consistency  of  the  present  model,  we  com¬ 
pare  in  the  next  section  the  predicted  saturated  and  nonsatur- 
ated  spectrum  with  measured  data. 


3.  INTEGRATION  TECHNIQUE  AND  SAMPLE  CALCULATIONS 


3 . 1  Method  of  Integration 

The  computer  code  that  has  been  developed  solves 
the  one  dimensional  wave  spectral  transfer  equation  (2.10.2) 
in  a  dimensionless  phase  space  whose  coordinates  are  C  and 
<.  The  code  uses  a  nonlinear  explicit  forward  marching  dif¬ 
ference  procedural  (in  time)  to  solve  the  equation.  These 
equations  are  written  in  a  coordinate  system  that  may  be 
moving  relative  to  the  ocean  and  in  which  it  is  assumed  that 
there  is  a  temporally  steady  but  spatially  variable  current 
[u(C)3.  The  velocity  at  which  energy  propagates  in  this 
coordinate  system  js  taken  to  be, 

U'  +  c  =  c^°Ma(k)  +  2/k(C)|  (3.1.1) 

relative  to  the  group  velocity  of  gravity  waves  in  a  static 
ocean . 

In  the  general  case  the  current  at  C  =  will  have 
a  finite  value  In  the  case  that  u_^  <  0,  the  energy 

propagation  velocity  will  be  positive  (to  the  right)  for 
sufficiently  small  wavenumbers  but  negative  for  very  large 
wavenumbers.  The  critical  wavenumber  (in  dimensionless  units) 
that  divides  these  two  regions  is 


(3.1.2) 


12  2 

k:  =  t  u  a  (<  )  for  u  <0  only 
c  4  oo  C  oo 

Boundary  conditions  when  <  0  roust  be  given  at  both 

^  =  +00  and  ^  =  -oo.  At  the  boundary  ^  =  -«>  values  for  the 

spectrum  Y  must  be  given  only  for  ;  whereas  at  C  = 

values  are  to  )^e  given  for  k>k  .  In  order  to  obtain  a  con- 

c 

vergent  solution  it  is  necessary  to  ensure  that  the  marching 
direction  is  always  in  the  direction  of  energy  propagation. 

At  any  station  C#  in  the  K,^-phase  space,  there  will 
be  a  critical  v;avenumber  which  separates  right  from  left 
propagating  waves, 

=  TU^(5)  '3-1- 

(see  Figure  6) .  Energy  is  transferred  across  this  boundary 
in  two  ways:  one,  by  local  gravity  wave-wave  interaction 
(modeled  by  the  diffusive  term  in  Eq.  (2.10.2)  and  two,  by 
convection.  Since  this  critical  boundary  corresponds  to 
zero  nu'c  group  velocity  in  physical  space,  the  trajectories 
of  wave  packets  in  the  plane  wi-1  have  a  vertical  slope 

at  this  boundary.  The  trajectories  (characteristics)  along 
which  energy  propagates  (in  the  absence  of  the  diffusive 
wave-wave  term)  are  given  (in  dimensionless  terms)  by. 
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It 


dx 


=  A(k) 
2/k 


+  u(U 


d< 

dx 


du  ) 


(3.1.4) 


The  family  of  trajectories  given  by  Eq.  (3.1.4)  for  the 
particular  current  pattern 


u(C) 


0.04 


evn  (-^ 


(3.1.5) 


is  depicted  in  Figure  (6b)  where  the  current  is  shown  in  (6a). 

In  the  absence  of  local  wave-wave  interactions  these 
curves  Fig.  (6b)  represent  the  direction  of  flow  of  energy 
in  the  (k,C) -phase  space.  A  wave  pacJcet  generated  at  some 
initial  point  will  travel  along  the  particular  trajectory 

which  passes  through  the  point  The  value  of  the  nor¬ 

malized  spoctrum  Y  along  this  trajectory  is  obtained  by  inte¬ 
grating  Eq.  (2.10.2)  in  the  form, 


dY 

37 


G(Y,k,C) 


(3.1.6) 


with 


dY  _  3Y  ^  3Y  dK  3Y 
37  3x  dx  T7  37  Tic 


(3.1.7) 


and 


37 


dK 

and  g-  are  the  trajectories  given  by  Eq.  (3.1.4) . 
Neglecting  the  diffusion  term  in  Eq.  (2.10.2)  has 
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essentially  decoupled  the  various  trajectories  except  for 
the  term  involving  A  from  Section  2.8.  Integration  along 
independent  characteristics  can  then  proceed  according  to 


As  mentioned,  to  ensure  stability,  the  direction  of  inte¬ 
gration  must  be  the  same  as  that  of  energy  flow.  Thus  the 
integration  starts  at  ^  for  and  follows  the 

energy  trajectories  in  Figure  (6);  while  for  <><  the 

crit 

integration  must  start  at  5  =  +oo. 


I 


'V 


3 . 2  Pierson  and  Stacy  Analytic  Fit  to  Data 

The  report  of  W.J.  Pierson  and  R.A.  Stacy  (1972) 
attempts  to  synthesize  the  experimental  data  on  the  develop¬ 
ment  of  surface  waves  as  a  function  of  wind  friction  velo¬ 
city  (U^).  One  major  difficulty  in  any  attempt  to  represent 
the  spectrum  of  the  ocean  surface  is  that  there  are  thirteen 
orders  of  magnitude  in  frequency  (six  in  wavenumber)  which 
must  he  de=:'ribed.  The  approach  taken  by  Pierson  and  Stacy 
(PS)  is  to  partition  the  wavenumber  spectrum  into  five  dis¬ 
tinct  regions  ani  to  fit  the  data  by  an  appropriate  spectral 
function  in  each  of  these  regions.  These  functions  are  then 
equated  at  the  boundaries  of  the  spectral  regions.  This 
phenomenological  approach  attempts  to  identify  the  major 
contributor  to  each  part  of  the  surface  spectrum  just  as  the 
spectral  transfer  equation  attempts  to  identify  the  inter¬ 
active  mechanism  in  the  equation  for  the  spectral  evolution 
(2.10.2).  The  empirical  expressions  of  PS  model  the  effects, 
whereas  Eq .  (2.10.2)  is  ar  attempt  to  model  the  causes. 

Comparison  of  calculations  with  the  empirical  expressions 
for  a  given  wind  speed  will  assist  us  Ln  determining  the 
veracity  of  our  models  of  the  individual  interactive  mechanisms. 

The  analytic  expression  constructed  to  represent  the 
spectrum  of  a  wind  roughened  sea  is  taken  directly  from  the  (PS) 
report.  We  first  define. 
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S(k,<I>)  =  S(k)  F(k,<I>)  0<k<-;-j<(t<^ 


(3.2.1) 


S(k,<I>)  =  0 


0^  k<'»;^<<ti<’!T  ;  -TT<il)<^ 


where  S(k)  is  the  one  dimensional  spectrum  discussed  in  this 
report  and  F(k,$)  cont-ains  all  the  angular  information, 


F(k,<l))  = 


8 


COS<I> 

3ti 


1  -  exp 


-g^/2k2[u(U*)] 


(3.2.2) 


,  1  ,  A  .  *  2. 

+  —  1  -  •-  +  A  cos  $ 

TT  2 


exp 


-  g^/2k^  [u(U*)  I 


where  0<A<2  and  A  could  be  a  function  of  k  and  U*  ,  S{k,^)  is 
the  wavenumber  spectrum  with  $  the  angle  between  )c  and  the  wind 
For  Tomparison  with  our  one  dimensional  theory,  i.e.,  the 
result  of  integrating  Eg.  (2.10.2),  we  set  F(k,<I>)  =  1. 

The  partitioning  of  the  wavenumber  specturm  is  done 
as  follows; 


S(k)  =  ^ 


(k) 

;  0  <  k  <  kj^ 

(i) 

1 

S2  (k) 

;  kj^  <  k  <  k2 

(ii) 

S3  (k) 

;  k2  <  k  <  k3 

(iii) 

(3,2.3) 

• 

S4(k) 

;  k,  <  k  <  k 

3  V 

(iv) 

I; 

d 

S3(k) 

8 

V 

yi 

V 

(V) 

The  analytic  forms  constructed  to  fit  the  experimental  data 
in  each  of  these  regions  are; 
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(i)  Gravity  wave-gravity  wave  equilibrium  range 


Si  (k)  =  — j  exp 


2k' 


2 


V/(Ui,_5(U.))  k 


(3.2.4a) 


0  <  k  <  k.  =  — = _ ni 


(ii)  isotropic  turbulence  range 


S_(k)  =  -  - 
2 


k^  <  k  <  k2 


(3.2.4b) 


(iii)  the  connecting  range  due  to  Leykin  and  Rosenberg 


a  D(U*) 
S,(k)  =  — — ^  - 
^  2kPk^"P 


k^  <  k  <  k^ 


(3.2.4c) 


(iv)  capillary  range 


S4(k)  = 


a  D(U*) 


2k- 


ks  <  k  <  k^ 


(3.2.4d) 


(v)  Viscous  cutoff  range 


S5(k)  = 


E  U*k® 
*  m 


vgk' 


k^  <  k  <  “> 


(3.2.4e) 


The  data  utilized  in  each  of  these  regions  has  not  been 
restricted  to  open  ocean  experiments;  use  has  also  been  made 
of  wind-water  tank  experiments.  The  bibliography  of  the 


experimental  literature  is  too  extensive  to  list  here;  see 
the  PS  report  for  this  information. 

In  addition  to  fitting  the  experimental  ’ata,  the 
functions  in  Eq.  (3.2.4)  were  selected  to  fit  smoothly  at 
the  boundciries  between  the  spectral  regions,  i.e.,  = 


S2(kj^)  , 


S2(k2)  =  S3(k2),  S3(k3)  =  S^Ck^)  and 


i\)  = 


The  parameters  in  Eq. 
following  values: 

a  =  8.1*10"^ 

6  =  0.74 
E  =  9.08  lO"^ 
E/vg  =  9.26  lO""^ 


(3.2.4)  were  determined  to  have  the 


k2  =  k  (f'Ti )  K  0.35 

k^  =  k(lOTT)  «  0.92 

k  =  (gp/T)*^  K  3.63 
m 

U*  =  12  cm/sec 
*m 


(3.2.5) 


and  functional  forms 

P  =  I 

and  (3.2.6) 

D(U*)  =  [l.274  +  0.0268  U*  +  6.03  10‘^U*p 


Ui9  ^(U*)  is  the  wind  speed  at  19.5  meters  for  a  friction 
velocity  of  U*  in  a  neutral  atmosphere.  A  complete  discussion 
relating  the  values  of  the  above  coefficients  to  the  appro¬ 
priate  experiments  may  be  found  in  the  (PS)  report. 

The  quantity  S(k)  defined  by  Eqs.  (3.2.^*)  corresponds 
to  the  equilibrium  value  of  the  one  dimensional  spectrum  F(k) 


calculated  from  the  spectral  transfer  equation  (2.10.2). 

In  Figures  (7)  and  (8)  we  graph  the  spectral  distribution  for 
two  values  of  the  wind  friction  velocity  U*  =  24  cm/sec  and 
48  cn./sec,  using  both  the  empirical  expressions  of  Pierson 

* 

and  Stacy  and  the  spectrum  obtained  by  integrating  Eq.  (2.10.2). 

3 

The  slope  spectrum  (Fk  )  has  been  graphed  vs.  the  log  of  the 
wavenumber  in  Figures  (7)  and  (8)  for  a  number  of  reasons: 
firstly,  because  the  spectra  span  decades  in  wavenumber; 
secondly,  because  the  Phillips'  saturated  spectrum  appears 
as  a  horizontal  line  in  this  representation;  thirdly,  the  sepa¬ 
ration  between  spectra  for  different  wind  riction  velocities 
is  most  clearly  made  and  finally,  because  of  the  clear  sepa¬ 
ration  between  the  regions  of  the  spectrum  as  modeled  by 
Pierson  and  Stacy,  it  is  an  easy  matter  to  compare  our  results 
with  theirs. 

In  Figure  (7)  we  see  that  the  low  wavenumber  (k  ^  .01) 
part  of  the  spectrum  is  modeled  equivalently  by  PS  or  the 
transfer  equation  for  a  U*  of  24  cm/sec.  The  gravity  wave- 
gravity  wave  equilibrium  range  |k^Sj^(k)j  essentially  rises 
to  the  value  a/2  where  the  PS  spectrum  L3comes  the  saturated 
spectrun  of  Phillips.  The  transfer  equation,  however,  deviates 
from  thfc  horizontal  giving  some  slope  to  this  region  -  the 
angle  of  which  is  dependent  on  U*.  The  saturation  is  there¬ 
fore  enhanced  above  the  value  given  by  PS.  The  transfer 
equation  rises  slowly  through  the  regions  given  by 

*  Values  of  the  free  parameters  0^,  E  and  were  chosen  to 
be  0.16,  0.5  and  0.0  for  this  comparison. 
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SLOPE  SPECTRUM  (Fk^) 


FIGURE  7. 


Fully-Developed  Spectrum  for  a  Friction 
Velocity  of  24  cm/sec. 


|k^S2(k)j  and  jk^S^Ck)  and  is  therefore  suppressed  below 
the  PS  spectrum  at  the  start  of  the  capillary  range.  If 
the  value  of  the  effective  viscosity  coefficient  (v^)  is 
taken  to  be  that  of  molecular  viscosity  the  transfer  spec¬ 
trum  becomes  too  large  in  the  capillary  regime.  This  growth 
represents  the  sensitive  balance  between  the  energy  input 
of  the  wind  and  energy  dissipation  due  to  viscous  damping 
in  Equation  (2.10.2).  With  modified  values  of  the  viscosity 
coefficient  =  (1. 7-2.5)  »  the  spectrum  is  bounded  in 

the  capillary  region.  It  is,  however,  clear  that  the  capil¬ 
lary  region  is  not  saturated. 


Figure  (8)  allows  us  to  make  two  separate  comparisons: 
first,  the  PS  spectra  for  two  values  of  the  wind  friction 
velocity,  i.e.,  U*  =  24  cm/sec  and  48  cm/sec  and  second,  how 
the  transfer  spectra  compare  with  the  PS  spectra.  We  see, 
first  of  all,  that  the  height  of  the  capillary  region  in  the 
PS  spectrum  scales  linearly  with  U*  ,  that  is,  for  twice 
the  wind  friction  velocity  one  has  twice  the  level  of  capillary 
waves.  In  fact  it  appears  that  outside  the  gravity  wave 
equilibrium  range  that  one  could  scale  the  PS  spectrum 
directly  with  U^.  Again,  for  the  transfer  spectrum,  (at 
U*  =  48  cm/sec ’/  the  balance  between  the  energy  input  of  the 
wind  and  the  energy  drain  of  viscosity  in  the  capillary 
|k^S^(k)j  and  cutoff  |k^S^(k)j  regions  produces  and  un- 
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bounded  growth.  This  growth  stresses  the  importance  of  con¬ 
structing  a  viable  model  for  the  viscosity  coefficient  which 
includes  a  dependence  on  the  wind  friction  velocity  U*.  An 
example  of  such  a  theory  might  be  the  application  of  the 
theory  of  breaking  waves  by  Banner  and  Phillips  (1973)  to 
determine  the  enhanced  viscosity  produced  by  turbulent  flow 
in  the  breaking  wave  region. 
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